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ABSTRACT 


A  compress ion- filter  design  is  described  which  utilizes  76  broadband 
taps  spaced  nonuniformly  along  a  delay  line.  Also  described  is  a  computer 
program,  written  in  FORTRAN  II  language,  which  can  be  used  to  calculate 
the  individual  tap  positions  and  weightings  for  a  maximum-response  peak- 
to-sidelobe  ratio.  This  program  is  quite  general  in  that  it  earn  be  used 
to  develop  an  optimum  (in  terms  of  response  peak-to-sidelobe  ratio)  com¬ 
pression  filter  for  an  arbitrary  filter  excitation-amplitude  function,  up 
to  100  taps,  and  any  time  delay-bandwidth  product  within  the  limitations 
of  the  computer  storage  capacity. 

Computations  made  using  the  computer  program  resulted  in  a  compres¬ 
sion  filter  matched  to  an  excitation  which  changed  frequency  linearly  with 
time  over  an  octave  bandwidth  and  which  had  a  time-bandwidth  product  of 
50.  The  response  of  this  filter  was  a  compressed  pulse  with  a  peak-to- 
sidelobe  ratio  of  over  100.  In  addition,  the  filter  proved  to  be  less 
sensitive  to  scan-rate  errors  than  a  continuously  dispersive  filter  with 
the  same  compression  characteristics. 

Experimental  results  were  obtained  on  a  76-tap  compression  filter 
designed  to  operate  in  the  1-  to  2-Gc  frequency  range.  The  performance 
of  this  filter  tends  to  confirm  the  validity  of  the  technique  and  provides 
an  encouraging  basis  for  future  work. 
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I .  INTRODUCTION 


A.  THE  PULSE- COMPRESS  ION  FILTER 

A  pulse-compression  filter  is  a  device  for  compressing  the  energy  of 
a  frequency-modulated  pulse  into  a  shorter  pulse  [Refs.  1  and  2],  Thus 
a  lossless  compression  filter  might  shorten  a  l-|isec,  1-v  (rms)  FM  pulse 
into  one  of  the  same  energy  with  a  duration  of  10  nsec  and  an  amplitude 
of  10  v.  The  degree  of  shortening  of  the  pulse  is  known  as  the  compres¬ 
sion  factor  (or  ratio)  and  would  be  100  for  the  above  filter.  The  com¬ 
pression  filter  operates  by  selectively  delaying  various  portions  of  the 
pulse  to  be  compressed  so  that  energy  at  the  beginning  of  the  pulse  is 
delayed  just  long  enough  to  emerge  from  the  filter  simultaneously  with 
energy  from  the  beginning  of  the  pulse. 

This  selective  delay  of  the  frequency-modulated  pulse  can  be  done  in 
several  ways.  One  method  is  to  use  a  network  whose  delay  varies  with 
frequency.  Another  is  to  use  bandpass  filters  to  separate  the  pulse 
energy  into  discrete  parts  according  to  frequency  and  then  to  recombine 
the  parts  through  the  use  of  suitable  delays  into  a  shorter  pulse  of 
larger  amplitude.  This  report  is  concerned  with  a  third  method--a  tapped 
delay  line  with  broadband  filters  on  each  tap. 


B.  THE  TAPPED  DELAY  LINE  WITH  BROADBAND  TAPS 


The  tapped-delay-line  compression  filter  is,  in  its  simplest  form,  a 
delay  line  with  provision  for  extracting  portions  of  the  signal  energy 
from  various  parts  of  the  line  and  adding  these  portions  together  at  the 
output.  This  is  illustrated  schematically  in  Fig.  1.  By  adjusting  the 
positions  of  the  taps,  the  various  portions  of  signal  energy  can  be  made 
to  add  in  phase  at  some  particular  time.  At  times  other  than  the  instant 
of  phase  addition,  some  degree  of  phase  cancellation  should  take  place 
depending  updn  the  number  of  taps,  their  relative  amplitude  contribution, 
and  the  characteristics  of  the  input  signal. 

Figur^  2  is  a  computer- generated  plot  of  the  response  of  a  tapped- 
delay-line  compression  filter  of  ten  taps  to  sinusoidal  excitation  of 


constant  magnitude  and  linearly  varying  frequency.  The  ten  sinusoids 
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representing  the  responses  of  the  individual  taps  are  shown  staggered  in 
time  by  amounts  equal  to  the  delays  of  the  taps.  The  bottom  trace  is 
the  sum  of  the  responses  of  the  ten  taps  and  is  seen  to  have  a  large  peak 
at  t  =  5  sec,  where  the  ten  taps  are  all  adding  in  phase.  The  other 
peaks  in  the  bottom  trace  are  reduced  in  amplitude  due  to  the  combined 
effects  of  fewer  taps  contributing  and  imperfect  phase  addition.  These 
secondary  peaks  are  analogous  to  grating  lobes  obtained  by  diffraction 
gratings  and  may  be  spaced  farther  apart  from  the  main  peak  by  increasing 
the  number  of  taps. 


C.  GENERAL  PROGRAM  DESCRIPTION 

The  computer  program  presented  in  this  report  was  designed  to  arrive 
at  the  optimum  delay  and  attenuation  of  each  tap  in  a  compression  filter 
to  obtain  the  maximum  ratio  of  the  amplitude  of  the  main  peak  to  the 
amplitude  of  the  signal  outside  the  region  of  the  main  peak.  This  was 
done  by  repeated  perturbations  of  the  individual  tap  positions  or  ampli¬ 
tudes  while  testing  for  improvement  in  peak-to-sidelobe  (PSL)  ratio  after 
each  perturbation.  The  sequence  was  arranged  to  operate  as  follows: 

1.  The  initial  tap  positions  and  amplitude  weightings  are  read  into 
the  program  as  data. 

2.  The  filter  response  is  calculated  with  the  initial  positions  and 
weightings . 

3.  The  amplitude  weighting  of  the  first  tap  (less  attenuation)  is 
increased . 

4.  The  filter  response  is  calculated  and  the  peak-to-sidelobe  ratio 
is  compared  with  the  original  value.  If  it  is  increased,  the 
operation  moves  to  the  next  tap.  If  it  is  not  increased,  the  sign 
of  the  perturbation  is  changed  and  the  test  is  repeated.  If  this 
change  results  in  improvement ,  the  tap  is  left  in  its  new  position 
and  the  next  tap  is  considered.  If  there  is  no  improvement  the 
tap  is  returned  to  its  original  position  before  moving  on. 

Z.  When  the  last  tap  is  reached,  the  final  peak-to-sidelobe  ratio  is 
compared  with  the  original  value.  If  the  improvement  is  less  than 
some  predetermined  value,  the  program  is  terminated.  Otherwise, 
the  above  sequence  is  repeated  using  delay  perturbations  instead 
of  weighting  perturbations. 


6.  The  sequences  of  amplitude  and  delay  perturbations  are  repeated 
alternately  with  the  peak-to-sidelobe  ratio  being  compared  each 
time  with  the  value  obtained  from  the  previous  sequence.  The 
program  terminates  when  the  increase  is  less  than  a  preset  value. 

The  program  was  written  in  FORTRAN  II  language  for  Stanford's  IBM 
7090  computer.  The  computational  portion  of  the  program  should  be  com¬ 
patible  with  other  computer  systems  comparable  to  the  7090  in  speed  and 
memory  capacity  (32,000  words).  Input  and  output  formats  may  have  to  be 
changed  to  run  this  program  on  other  systems;  the  plot  routines  will 
almost  undoubtedly  have  to  be  changed  since  those  used  are  unique  to 
the  Stanford  computer. 

The  flow  chart  of  Appendix  A  should  provide  a  useful  device  for  the 
reader  desiring  to  obtain  a  knowledge  of  the  mechanics  of  the  computer 
program.  It  should  be  especially  helpful  to  a  programmer  in  changing 
the  program  language  or  in  making  any  other  program  modifications. 

D.  MATHEMATICAL  DESCRIPTION  OF  PROGRAM 

In  each  of  the  computer  runs  the  filter  excitation  is  expressed  by 
the  following  equation: 

e(t)  =  A(t)  cos  [2n(2t  -  O.Olt2)]  (l) 

This  equation  represents  a  signal  whose  frequency  is  changing  as  a  linear 

function  of  time.  The  normalized  times  and  frequencies  used  in  this 

program  were  chosen  for  computational  convenience  and  could,  for  instance, 

represent  nanoseconds  and  gigacycles  as  well  as  seconds  and  cycles. 

th 

The  response  of  the  n  tap  to  the  above  excitation  is 

rn(t)  =  A(t-Tn)  Bn  cos  ^2jx  j^2( t-Xn )  -  0.0l(t-Tn)2j|  (2) 

where 

=  time  delay  to  the  nt*1  tap  from  the  filter  input  (sec) 

B  =  amplitude  multiplying  factor  of  the  n**1  tap 


SEL-65-043 


4 


The  response  of  a  complete  filter  of  m  taps  is  the  sum  of  the  responses 
of  the  individual  taps: 


f 


m 

r(t)  =  ^  rn(t)  (3) 

n=l 


The  coefficients  in  the  above  equations  were  chosen  to  yield  a  time- 
bandwidth  product  of  50  over  an  octave  bandwidth,  as  may  be  seen  by 
taking  the  time  derivative  of  the  phase  of  Eq.  (l)  to  obtain  the  fre¬ 
quency  : 


I 


I 


{ 


I 


oj(t)  =  ^  [«(t)]  =  2n(2  -  0.02t)  (4) 

The  frequency  is  2  cps  at  t  =  0  sec  and  will  be  1  cps  at  t  =  50  sec. 

If  the  values  of  Tq  are  chosen  such  that  the  outputs  of  all  the  taps 
add  in  phase  at  t  =  50  sec,  then  the  compressed  pulse  will  have  the 
desired  time  delay-bandwidth  (time-bandwidth)  product  of  50.  The  proper 
values  of  Tq  are  obtained  by  setting  t  s  50  sec  in  Eq.  (2)  and  solving 
for  the  values  of  which  will  cause  the  cosine  argument  to  equal  inte¬ 

gral  multiples  of  2n: 

2n|^2(50-Tn)  -  0.0l(50-Tn)2]  =  2nk  (5) 

where  k  is  an  integer.  Solving  for  T^, 

Tk  =  -50  +  /10,000  -  100k  (6) 

Equation  (6)  indicates  that  there  are  76  values  of  k  from  0  to  75  that 
can  be  used  to  calculate  usable  T's  from  0  to  50.  If  n  is  taken  to 
be  the  tap  number  counting  from  the  filter  input,  then  k  =  75-n  and 
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T  =  -50  +  J 2500  +  lOOn  (7) 

n  v 

These  are  the  initial  values  of  which  are  used  in  Eq.  (3)  and  about 

which  the  perturbations  are  made. 

Provision  is  made  in  the  program  to  set  the  initial  values  of  B 

n 

at  any  desired  level.  The  magnitudes  of  the  perturbations  may  be  chosen 
at  any  initial  values  and,  if  desired,  the  perturbations  may  be  changed 
to  smaller  values  in  midprogram.  This  might  be  desirable  since  it  allows 
the  program  to  converge  on  a  solution  with  fewer  steps. 

The  function  A(t)  in  Eq.  (l)  represents  the  composite  bandpass 
magnitude  function  of  all  elements  in  the  rf  signal  path,  whether  the 
elements  precede  or  follow  the  filter.  In  a  microwave  filter  of  octave 
bandwidth,  the  magnitude  function  will  probably  be  determined  predomi¬ 
nately  by  external  components ,  such  as  an  associated  TWT  amplifier. 

Several  functions  are  available  in  the  program  in  the  form  of  subroutines, 
but  the  main  optimizing  program  used  a  parabolic  shaping  function  which 
approximated  the  passband  shape  of  the  WT  that  was  used  to  test  the 
prototype  filter. 

In  all  portions  of  the  program,  time  was  incremented  in  steps  of  0.1 
sec.  This  value  was  chosen  as  a  compromise  between  computational  effi¬ 
ciency  and  accuracy.  The  average  amplitude  error  introduced  by  this 
increment  size  is  about  3  percent  in  the  sidelobes  where  the  time  samples 
are  random  relative  to  the  peaks,  and  about  1  percent  in  the  compressed 
pulse  which  is  initially  fixed  relative  to  the  time  samples.  Since 
improvements  in  peak-to-sidelobe  ratios  of  several  hundred  percent  were 
sought,  it  was  not  expected  that  these  errors  would  have  any  significant 
influence  on  the  results. 

The  calculation  of  a  PSL  ratio  requires  that  the  computer  be  able  to 
distinguish  between  a  peak  and  a  sidelobe.  Normally,  the  first  envelope 
nulls  on  either  side  of  the  main  peak  might  be  chosen  as  the  boundary 
between  "peak"  and  "sidelobe."  The  largest  peak  between  nulls  would  then 
be  used  to  determine  the  PSL  ratio.  This  technique  was  quickly  deter¬ 
mined  to  be  impractical  when  preliminary  filter-response  plots  obtained 
on  the  computer  showed  that  clearly  defined  envelope  nulls  did  not 
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necessarily  exist  in  the  region  of  the  main  peak.  In  addition,  the 
simulation  of  a  detector  increased  computer  running  time  considerably. 

An  alternative  method  of  choosing  the  region  of  the  main  peak  is 
possible  due  to  the  fact  that  the  approximate  location  in  time  of  this 
peak  is  known  from  Eq.  (5).  This  alternative  simply  involves  having  the 
computer  look  for  the  sidelobes  outside  some  fixed  time  interval  sur¬ 
rounding  the  time  chosen  in  Eq.  (5)  in  order  for  phase  addition  to  occur. 
The  optimization  program  described  here  used  a  time  interval  of  ±3  sec. 
Inspection  of  the  computer  plots  of  Chapter  IV  will  show  that  sidelobe 
suppression  outside  this  time  interval  was  not  obtained  at  the  expense 
of  larger  peaks  immediately  inside  the  ±3  sec  boundaries.  It  appears 
that  the  program  will  be  insensitive  to  the  exact  locations  of  these 
boundaries  if  they  do  not  fall  within  the  body  of  the  main  pulse. 

E.  INFORMATION-OUTPUT  TECHNIQUE 

The  availability  of  the  tape-fed  plotter  made  the  amplitude- time 
plot  a  natural  form  of  output  for  this  program.  Not  only  does  the  plot 
reduce  the  bulk  and  ease  the  interpretation  of  the  output,  but  it  pro¬ 
vides  the  same  display  that  one  would  observe  on  an  oscilloscope  of  an 
actual  filter  output. 

All  plots  were  scaled  to  10  sec  per  inch  on  the  horizontal  scale. 

The  vertical  scale  was  chosen  such  that  the  absolute  magnitude  of  the 
compressed  pulse  was  4.5  in.,  this  value  being  chosen  so  that  the  plot 
would  fit  onto  a  standard  notebook  page.  An  expanded  plot  with  the 
vertical  amplitude  scale  multiplied  by  10  was  provided  to  allow  closer 
examination  of  the  sidelobe  structure.  Calibration  tick  marks  were 
provided  at  the  side  of  the  plot  at  0.1,  0.05,  0.025,  0.0125,  and  0.01 
times  the  peak  amplitude. 

The  printed  output  includes  the  location  of  each  T  ,  the  magnitude 
of  each  B^,  the  magnitude  of  the  main  peak  and  its  location,  the  mag¬ 
nitude  of  the  largest  sidelobe  and  its  location,  the  PSL  ratio,  and  the 
improvement  in  PSL  ratio  obtained  by  the  last  series  of  perturbations. 
This  information  is  printed  each  time  the  computer  finishes  the  series 
of  perturbations  of  the  m  taps,  thus  providing  an  indication  of  the 
rate  at  which  improvement  is  being  made  in  the  filter. 
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II.  COMPUTATIONAL  RESULTS 


A.  DESIGN  OBJECTIVES 

The  computer  program  was  applied  to  the  design  of  a  microwave  com¬ 
pression  filter  for  use  with  high  scan-rate  excitation.  The  desired 
parameters  were:  bandwidth  =  1  Gc,  center  frequency  =  1.5  Gc,  and  time- 
bandwidth  product  =  50.  The  desired  weighting  function  was  approximated 
by  the  following  parabolic  function: 

for  -10  g  t  g  60 

(8) 

l  0  for  -10  >  t  >  60 

where  t  is  in  nanoseconds.  Then  the  desired  excitation  function  is 

10  g  t  g  60 

(9) 

10  >  t  >  60 

B.  OPTIMIZATION  PROGRAM  OUTPUT  DATA 
1 .  Computer  Run  A 

a.  Initial  Conditions 
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B  =  0.5 
n 


.  m  =  76  taps;  n.=  0,  1,  2,  ....  74,  75 

,  T  =  -50  +  100^1.01  -  (0.0l)n 

n 

AB  =  0.1  (perturbation  increment) 

At  =  0.1  (perturbation  increment) 

Cutoff  improvement  =  1  percent  of  PSL  ratio 

Time  cutoff  =  30  min  (computer  allows  this 
value  to  be  exceeded  by  10 
percent  before  terminating) 

b.  Output 

Run  A  terminated  after  being  unable  to  obtain  any  improvement 
in  PSL  ratio  on  the  fourth  series  of  perturbations  of  the  T's.  The 
improvements  obtained  on  each  series  are  shown  in  Fig.  3.  Total  running 
time  of  the  program  on  the  7090  computer  was  24.0  min,  of  which  3.6  min 
was  compile  time. 

Remarks :  Almost  all  of  the  improvement  in  this  filter  was 
obtained  through  perturbation  of  the  B's.  In  addition,  the  overall 
trend  of  the  perturbations  was  to  increase  the  values  of  the  B's  in 
the  center  of  the  filter  while  decreasing  them  at  the  ends.  It  was 
decided  to  rerun  the  program  perturbing  only  the  B's  and  setting  an 
initial  weighting  on  their  magnitudes. 

I 

2.  Computer  Run  B 

[ 

i 

a.  Initial  Conditions 

The  initial  conditions  for  this  program  were  the  same  as  for 
•  run  A  except  that  At  was  not  used,  AB  =  0.05,  and  the  B's  were  as 

shown  in  Table  1. 
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PERTURBATION  SERIES 

BrBrBrBr 
VARIABLE  PERTURBED 

34910 


FIG.  3.  PEAK-TO-SIDELOBE 
RATIO  AND  IMPROVEMENT  IN 
PSL  RATIO  FOR  EACH  SERIES 
OF  PERTURBATIONS  OF  THE 
76  TAPS,  COMPUTER  RUN  A. 


TABLE  1.  INITIAL  VALUES  OF  B  FOR  COMPUTER  RUN  B 

n 


Tap 

0-2 

3-6 

7-11 

17-20 

21-25 

25-30 

31-34 

35-40 

Number 

73-75 

69-72 

64-68 

55-58 

50-54 

45-49 

41-44 

B 

n 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

b.  Program  Modifications 

The  section  of  Program  II  (see  Appendix  B)  which  is  respon¬ 
sible  for  the  perturbations  of  the  T's  was  removed.  This  is  the  portion 
of  the  program  between  the  arrows  on  th£  program  listing. 

c.  Output 

This  run  terminated  after  obtaining  only  a  0. 46-percent 
improvement  in  PSL  ratio  on  the  fifth  perturbation  series.  The  fluctua¬ 
tion  of  the  improvements  (Fig.  4)  indicated  that  more  running  time  with 
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a  smaller  cutoff  change  percentage  might  be  fruitful.  Running  time  was 
13  min. 


[—RUN  B— |  1- - —RUN  C 


PERTURBATION  SERIES 

34912 

FIG.  4.  PSL  RATIO  AND  IMPROVEMENT  FOR 
COMPUTER  RUNS  B  AND  C. 

3.  Computer  Run  C 

a.  Initial  Conditions 

The  cutoff  change  percentage  was^  changed  to  0.2  percent  and 
a  limit  of  12  was  put  on  the  number  of  perturbation  series  allowed.  The 
perturbation  increment  was  changed  to  0.01,  and  the  initial  B's  were 
those  arrived  at  in  rim  B. 

b.  Output 

The  program  terminated  after  12  perturbation  series  were 
performed  without  encountering  a  change  in  PSL  ratio  less  than  0  2  per¬ 
cent.  The  ratios  and  improvement  percentages  are  plotted  in  Klg.  4,  and 
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the  migrations  of  the  weightings  of  the  76  taps  are  plotted  in  Fig.  5 
for  both  runs  B  and  C.  Figures  6  and  7  are  plots  of  r(t)  vs  time  for 
the  final  set  of  B's  and  T's,  whose  values  are  given  in  Table  2. 

The  running  time  was  25  min. 


C.  COMPUTER  TESTS  OF  OPTIMIZED  FILTER 


Sensitivity  of  Filter  to  Random  Excitation  Phase 


Program  I  was  used  to  calculate  the  response  of  the  filter 
derived  in  run  C  to  an  arbitrary  phase  shift  in  the  input  function.  The 
ten  phase  values  used  were  obtained  from  a  random  number  table.*  The 
function  calculated  is 


76 

r  (t)  =  Y  A(t-T  )  Bn  cos  2n|2(t-T  )  -  O.Ol(t-T  )2  +  ♦  1  (ll) 
i  Z_,  n  n  L  n  n  ij 

n=l 


where 

i 


-10  §  (t-T  )  §  60 
n 


-10  >  (t-T  )  >  60 
n 


B  and  T  are  as  given  in  Table  2,  and 
n  n 

*1  =  0.10480  04  =  0.02011  U>7  =  0.69179 

=  0.15011  $_  =  0.81647  ="0.  J.4194 

/do 

=  0.01536  =  0.91646  ♦_  =  0.62590 

0  0  7 

*10  =  °- 36207 

* 

C.R.C.  Standard  Mathematical  Tables,  Chemical  Rubber  Company,  Cleveland, 
Ohio,  1959. 
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FIG.  6.  OPTIMIZED  FILTER  RESPONSE  (COMPUTED). 


0  10  20  30  40  30  60  70  80  90  100 

TIME  (sec) 

FIG.  7.  OPTIMIZED  FILTER  RESPONSE  (VERTICAL  SCALE  EXPANDED  10  TIMES). 

The  ten  plots  obtained  were  traced  over  each  other  (by  hand),  which 
resulted  in  the  composite  plot  of  Fig.  8. 

2.  Filter  Response  for  Various  Scan  Rates 

To  test  the  sensitivity  of  the  optimized  filter  to  deviations 
from  the  design  scan  rate,  the  coefficient  of  the  quadratic  time  term 
was  changed  in  Subprogram  ESUMI  and  this  revised  function  was  then  used 
with  Program  I  to  calculate  the  filter  response.  The  following  coef¬ 
ficients  were  used:  0.005,  0.00P,  0.0095,  0.01,  0.0105,  0.011,  0.012, 
0.015,  0.02.  These  correspond  to  scan  rates  of  10,  16,  18,  19,  20,  21, 
22,  24,  30,  and  40  cycles/sec2  respectively.  The  results  of  these  tests 
are  plotted  in  Fig.  9. 

3 .  Effect  of  Various  Excitation  Function  Envelopes 

The  four  functions  of  Fig.  10  were  chosen  to  test  the  effect  on 
the  optimized  filter  of  a  wide  variation  in  input  function  shapes.  The 
four  functions  used  are  shown  beginning  on  page  20. 
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TABLE  2.  VALUES  OF  B  AND  T  FOR  76-TAP  FILTER  WITH  PSL  RATIO  =  107 


Tap 

Number 

B 

T 

(nsec) 

1 

0.08 

0.000 

2 

.12 

0.990 

3 

.04 

1.962 

4 

.10 

2.915 

5 

.14 

3.852 

6 

.14 

4.772 

7 

.14 

5.678 

8 

.19 

6.568 

9 

.27 

7.446 

10 

.23 

8.310 

11 

.29 

9.161 

12 

.28 

10.000 

13 

.30 

10 . 828 

14 

.31 

11.644 

15 

.31 

12.450 

16 

.35 

13.246 

17 

.31 

14.031 

18 

.33 

14.807 

19 

.35 

15.574 

20 

.33 

16.332 

21 

.36 

17.082 

22 

.36 

17.823 

23 

.33 

18.556 

24 

.37 

19.282 

25 

.41 

20.000 

Tap 

Number 

B 

T 

(nsec) 

Tap 

Number 

B 

26 

0.37 

20 

.711 

51 

0.45 

27 

.41 

21 

.414 

52 

.39 

28 

.43 

22 

.111 

53 

.39 

29 

.44 

22 

.801 

54 

.35 

30 

.41 

23 

.485 

55 

.36 

31 

.42 

24 

.162 

56 

.31 

32 

.46 

24 

.833 

57 

.30 

33 

.46 

25 

.498 

58 

.33 

34 

.46 

26 

.158 

59 

.29 

35 

.52 

26 

.811 

60 

.26 

36 

.47 

27 

.460 

61 

.25 

37 

.52 

28 

.102 

62 

.21 

38 

.45 

28 

.740 

63 

.21 

39 

.51 

29 

.372 

64 

.24 

40 

.51 

30 

000 

65 

.15 

41 

.52 

30 

.623 

66 

.15 

42 

.49 

31 

240 

67 

.14 

43 

.51 

31 

853 

68 

.15 

44 

.48 

32 

462 

69 

.15 

45 

.47 

33 

066 

70 

.10 

46 

.50 

33 

666 

71 

.11 

47 

.42 

34 

261 

72 

.10 

48 

.51 

34 

852 

73 

.10 

49 

.42 

35 

440 

74 

.10 

50 

.43 

36 

023 

75 

.05 

T 

(nsec) 

36.602 
37.178 
37.750 
38.318 
28.882 
39.443 
40.000 
40.554 
41.104 
41.651 
42.195 
42.736 
43.274 
43.808 
44.340 
44 . 868 
45.394 
45.917 
46.436 
46.954 
47.468 
47.980 
48.489 
48.995 
49.499 
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50 . 000 


20  30  40  50  60  70  80  90 

TIME  (toe) 

FIG.  8.  COMPOSITE  PLOT  OF  RESPONSE  OF  OPTIMIZED  FILTER  TO 
EXCITATION  WITH  RANDOM  PHASE  CONSTANT. 
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FIG.  10.  FOUR  FUNCTIONS  USED  AS 
EXCITATION  AMPLITUDE- WEIGHTING 
FUNCTION,  A(t). 
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1 


for 


-10  s  (t-x)  ^  60 


( 


t-T-25 \2 
35  / 


A ( t-x)  =  { 


(12) 


for  -10  >  (t-x)  >  60 


A (t-X )  =  i 


|cos  n  | 


^)f 


for  -10  §  (t-x)  §  60 


(13) 


for  -10  >  (t-x)  >  60 


A(t-X)  = 


£cos  It  ^ 


^)]1/2 


for  -10  §  (t-x)  §  60 


(14) 


for  -10  >  (t-x)  >  60 


A(t-x)  =  { 


[cos  „  (i^)] 


for  -10  §  (t-x)  g  60 


(15) 


for  -10  >  (t-x)  >  60 


The  compressed-pulse  parameters  resulting  from  these  four  functions  are 
listed  in  Table  3. 
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TABLE  3.  FILTER  PERFORMANCE  FACTORS 
FOR  DIFFERENT  EXCITATION  MAGNITUDE  FUNCTIONS 


Function 

6-db 

Pulse  Width 
(sec) 

PSL 

Ratio 

Compressed -Pulse 
Amplitude 

t2 

1 - 2 

1.4 

107 

21.2 

352 

(cos  x)2 

1.7 

100 

18.5 

(cos  x)1/2 

1.3 

61 

22.1 

(cos  x)1/4 

1.3 

40 

22.9 

4 .  Image  Response 

The  sign  of  the  scan  rate  was  changed  in  Subprogram  ESUMI  to 
determine  the  image  response  of  the  optimized  filter.  The  image  excita¬ 
tion  function  is: 


e(t)  =  A ( t )  cos  2n(l.0t  +  O.Olt2)  (16) 

where  A(t)  is  given  by  Eq.  (8).  The  calculated  response  is  plotted  in 
Fig.  11  to  the  same  time  and  amplitude  scale  as  the  optimized  filter 
response  of  Fig.  6.  The  peak  magnitude  of  the  image  was  14.8  db  below 
that  of  the  main  pulse. 

D.  REMARKS  ON  COMPUTATIONAL  RESULTS 

The  early  termination  of  run  A  was  due  to  the  inability  of  the  pro¬ 
gram  to  improve  the  peak- to-sidelobe  ratio  by  perturbing  the  T's.  As 
can  be  seen  in  Fig.  3,  none  of  the  four  series  of  T  perturbations 
resulted  in  any  significant  improvement  in  the  filter.  Since  the  migra¬ 
tions  of  the  B's  in  run  A  indicated  a  trend  toward  something  resembling 
a  Taylor  weighting  [Ref.  l],  a  stepped  first-order  approximation  to  this 
weighting  was  used  in  run  B  for  the  initial  B's  in  order  to  decrease 
the  computation  time. 
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FIG.  11.  RESPONSE  OF  OPTIMIZED  FILTER  TO  IMAGE  EXCITATION. 

The  fluctuations  of  the  weightings  as  shown  in  Fig.  5  indicate  that 
the  one-at-a-time  perturbation  technique  is  not  the  most  efficient  method 
of  deriving  an  optimum  filter.  It  does,  however,  have  the  advantage  of 
simplicity.  If  a  new  program  were  to  be  written,  it  would  probably  be 
better  to  perturb  the  taps  in  variable-size  groups  instead  of  individually. 

The  finally  derived  filter,  with  a  PSL  ratio  of  107,  allowed  a  50- 
percent  scan-rate  error  before  its  PSL  ratio  was  degraded  to  that  of  the 
constant-weighting  filter  (Fig.  3).  The  pulse  width  is  increased  by 
about  a  factor  of  10  with  a  50-percent  scan-rate  error  due  to  the  much 
smaller  bandwidth  over  which  phase  addition  can  take  place. 
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III.  EXPERIMENTAL  RESULTS 


As  a  parallel  Investigation  to  the  computer  synthesis  of  a  irpped- 
delay-line  compression  filter,  the  realizability  of  such  a  fi-iter  at 
microwave  frequencies  was  explored.  The  purpose  of  constructing  this 
filter  was  to  determine  if  phase  and  attenuation  could  be  controlled 
to  the  required  accuracy  in  a  microwave  device  with  as  many  as  76  taps. 

In  accordance  with  the  numbers  used  in  the  optimization  program,  the 
filter  was  designed  to  operate  between  1  and  2  Gc,  with  a  time-bandwidth 
product  of  50. 

A .  CONSTRUCTION 

The  delay  line  was  constructed  from  shielded  stripline  due  to  the 
compactness  of  the  structure  and  the  accessibility  of  the  conductor  for 
coupling  purposes.  It  consists  of  two  parallel  strips,  each  237  in.  - 
long,  photoetched  onto  a  Tellite  3B  dielectric  sheet  0.0625  in.  thick. 
Another  dielectric  sheet  of  the  same  thickness  is  laid  on  top  of  the 
conductors,  and  the  two  dielectric  sheets  are  then  sandwiched  between 
two  aluminum  ground  planes  to  form  the  complete  delay  structure.  The 
delay  line  with  one  ground  plane  and  dielectric  sheet  removed  is  shown 
in  Fig.  12.  The  large  number  of  screw  holes  was  necessary  to  maintain 
the  spurious  coupling  level  between  the  parallel  lines  at  a  value  lower 
than  -50  db.  A  drawing  of  the  line  is  shown  in  Fig.  13. 

The  couplers  can  be  resistive  or  reactive  devices  connecting  the  two 
parallel  lines  at  the  appropriate  points  to  achie/e  the  delays  given  in 
Table  2.  Since  signals  on  the  output  line  can  ccuple  back  into  the  input 
line  and  thence  into  the  output  line  again  through  any  combination  of 
couplers,  it  is  necessary  that  the  couplers  have  at  least  20-db  loss  to 
prevent  spurious  signals  from  appearing  at  the  output  in  the  form  of 
sidelobes.  This  requirement  of  high  coupler  loss  could  be  relaxed  some¬ 
what  if  directional  couplers  were  used,  but  with  76  couplers  the  limit 
is  still  about  -15  db  if  sufficient  power  is  to  be  conserved  for  the  last 
tap . 
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FIG.  12.  PHOTOGRAPH  OF  PARALLEL  TRANSMISSION  LINES 
BONDED  TO  TELLITE  3B  SUBSTRATE. 

It  was  initially  planned  to  use  resistive  taps  between  the  two  sec¬ 
tions  of  the  delay  line.  However,  the  fact  that  the  resistors  necessarily 
remove  considerably  more  power  from  the  input  line  than  is  delivered  to 
the  output  line,  makes  this  type  of  coupler  unusable  in  this  application 
because  of  the  resulting  high  attenuation.  The  use  of  edge-coupled 
quarter-wavelength  directional  couplers  [Ref.  3]  was  a  logical  considera¬ 
tion,  but  the  construction  and  placement  of  76  such  couplers  would 
involve  considerable  expense  and  time.  The  siqpler  coupler  configura¬ 
tion  actually  used  was  intended  primarily  to  demonstrate  whether  phase 
could  be  accurately  controlled  through  76  different  signal  paths  simul¬ 
taneously  . 

Figure  14  is  a  drawing  of  the  coupler  configuration  used,  and  Fig. 

15  is  a  plot  of  the  frequency  response  of  this  coupler.  It  is  simply  a 
piece  of  No.  40  wire  laid  across  the  conductors  perpendicularly  and 
separated  from  them  by  a  0.002-in,  piece  of  cellophane  tape.  Strips  of 
0.007-in.  cellophane  tape  were  laid  between  the  noncoupled  conductors  to 
serve  as  spacers  and  to  prevent  the  copper  wire  from  It ing  squeezed 
through  the  0.002- in.  tape. 
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COPPER  CONDUCTORS 
0.002  X  0.093 


34905 


FIG.  14.  CUTAWAY  SIDE  VIEW  OF  COUPLER  CONFIGURATION. 


l.O  1.9  2.0 


FREQUENCY  (Gc) 

34909 

FIG.  15.  FREQUENCY  RESPONSE  OF  THE 
COUPLER  CONFIGURATION  OF  FIG.  14. 
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The  locations  of  the  taps  along  the  delay  line  were  determined  by 
the  values  of  the  T's  in  Table  2.  Tap  1  was  used  as  a  reference  with 
T  =  0,  and  the  total  length  of  delay  line  from  tap  1  to  tap  n  is 
given  by: 


1  ,n 


inches 


(17) 


where  e  is  the  relative  dielectric  constant  of  the  substrate,  (e  = 
r  r 

2.32  for  Tellite  3B),  and  11.803  in. /nsec  is  the  speed  of  light  in  vacuo. 

The  distance  d,  doesn't  include  the  length  of  the  coupler  since  that 
1  ,n 

is  a  constant  value  for  each  tap.  The  electrical  length  of  the  corners 
was  determined  empirically  to  be  0.0147  nsec  or  0.114  in.  for  each  corner. 
The  dimensions  required  for  a  given  delay  are  clarified  in  Fig.  16. 


B.  PERFORMANCE 

The  test  apparatus  of  Fig.  17  was  used  to  examine  the  compression 
characteristics  of  the  prototype  filter.  Photographs  of  the  filter 
response  are  shown  in  Figs.  18  through  22.  These  photographs  do  not 
necessarily  represent  the  full  capabilities  of  the  filter  because  of 
spurious  responses  from  the  TWT's  and  scan-rate  variations  in  the  BW0. 

The  dual-trace  photograph  of  Fig.  18  shows  the  compressed  and  uncom¬ 
pressed  rf  envelopes  in  their  proper  time  relationships.  The  uncompressed 
pulse  at  the  filter  input  is  attenuated  by  20  db  in  the  photograph,  thus 
indicating  a  net  compression  loss  of  20  db.  From  Fig.  15  the  filter  loss 
at  the  WT  center  frequency  (1300  Me)  is  34  db,  which  gives  a  compression 
ratio  of  approximately  14  db  or  25  to  1 .  This  result  suggests  that  the 
phasing  of  the  individual  taps  is  accurate  since  an  in-phase  contribution 
from  almost  all  of  the  taps  is  required  to  achieve  a  14-db*  compression 
ratio  with  the  limited  bandwidth  of  the  TWT's.  The  magnitude  of  the 

* 

The  compression  factor  is  equal  to  the  product  of  the  bandwidth  over 
which  in-phase  contributions  are  obtained  times  the  differential  time 
delay  of  that  bandwidth  in  the  compression  filter.  This  "time-bandwidth 
product"  is  explained  in  detail  in  Ref.  1. 
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34907 

FIG.  16.  DIAGRAM  SHOWING  METHOD  OF  MEASURING  ELECTRICAL  DISTANCE  BETWEEN 

TWO  TAPS.  In  general,  d  =  (d  +  d  +  d  +  d  +  d  +  d  )  +  4d  . 

K  |  K*  1  x  a  J  ■  *)  b  ( 

For  the  dimensions  shown,  d?  =  0.114  in.,  where  d?  is  the  electrical 
distance  around  a  corner. 


34913 

FIG.  17.  SIGNAL  FLOW  DIAGRAM  OF  APPARATUS  FOR  TESTING  THE  76-TAP  FILTER. 
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FIG.  18.  DUAL-TRACE  PHOTOGFAPH  OF  76-TAP 
FILTER  RESPONSE  AND  EXCITATION.  Lower 
trace:  compressed-pulse  envelope,  10 

nsec/div.  Upper  trace:  excitation 
envelope,  same  relative  time  scale. 

Amplitude  is  proportional  to  power  in 
both  traces. 

frequency  response  of  the  TWT's  is  the  shape  of  the  uncompressed  pulse 
of  Fig.  18,  with  frequency  decreasing  from  left  to  right  and  the  point 
directly  over  the  compressed  pulse  corresponding  to  1  Gc.  The  10-nsec/div 
time  scale  converts  to  a  200-Mc/div  frequency  scale  as  a  result  of  the 
20-Mc/nsec  scan  rate. 

The  "spreading"  effect  of  thte  compression  filter  upon  a  signal  scan¬ 
ning  in  the  wrong  direction  (such  as  an  image)  is  shown  in  Fig.  19.  The 
image  energy  is  spread  out  in  time  to  approximately  twice  its  original 
duration,  giving  the  signal-to- image  ratio  as  the  compression  ratio  times 
2,  or  in  this  case,  14  +  3  =  17  db.  This  ratio  is  not  apparent  in  Fig. 

19  due  to  the  fact  that  the  sampling  oscilloscope  distorts  signals  whose 
durations  are  very  short  compared  to  the  time  scale. 

The  duration  of  the  compressed  pulse,  when  observed  on  the  expanded 
time  scale  of  Fig.  20,  is  seen  to  approach  that  of  the  impulse  response 
of  the  WT's  (Fig.  21 ).  For  comparison  purposes,  it  should  be  noted  that 
the  compressed  pulse  is  a  detected  envelope,  while  the  actual  rf  cycles 
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of  the  T\VT  impulse  response  are  shown.  Also,  the  compressed-pulse  ampli¬ 
tude  is  proportional  to  power,  while  the  impulse-response  amplitude  is 
proportional  to  voltage. 


FIG.  19.  COMPRESSED  PULSE  AND  IMAGE  RESPONSE, 
25  NSEC/DIV.  Amplitude  is  proportional  to 
power. 


FIG.  20.  COMPRESSED- PULSE  ENVELOPE ,  2  NSEC/DIV. 
Amplitude  is  proportional  to  power. 


FIG.  21.  IMPULSE  RESPONSE  OF  TWT  AMPLIFIERS 
USED  IN  COMPRESSION  SYSTEM,  2  NSEC/DIV. 


The  last  photograph,  Fig.  22,  illustrates  the  sidelobes  present  along 
the  baseline  of  the  compressed  pulse.  These  sidelobes  are  caused  by  a 
combination  of  the  following  effects: 

1.  Output  power  variations  in  the  BWO  contain  significant  L-band 

components  due  to  the  rapid  scan  rate.  The  mixer  converts  these 
variations  into  energy  within  the  compression  filter  passband , 
where  it  appears  as  baseline  clutter. 


FIG.  22.  COMPRESSED  PULSE  AND  SIDELOBES,  10 
NSEC/DIV.  Main  pulse  amplitude  =  10  divisions. 


2.  The  TWT's  used  in  this  test  have  internal  reflections  which  cause 
delayed  and  attenuated  replicas  of  an  input  pulse  to  appear  at 
the  output  for  up  to  100  nsec  following  the  direct  transmission. 

3.  Imperfect  phase  cancellation  of  the  compression  filter  results  in 
peaks  of  energy  outside  the  region  of  the  compressed  pulse. 

4.  The  instantaneous  frequency  of  the  scanning  excitation  deviates 
from  linear  variation  by  ±10  Me  at  several  points  in  the  sweep. 

Most  of  the  sidelobes  of  Fig.  22  are  caused  by  effects  3  and  4.  The 
two  largest  ones,  however,  are  due  to  effects  1  and  2.  The  peak  preceding 
the  pulse  by  50  nsec  is  due  to  direct  detection  and  is  20  db  lower  than 
the  main  pulse.  The  peak  following  the  main  pulse  by  20  nsec  is  also  down 
20  db  and  is  due  to  TWT  reflections.  The  sidelobes  due  to  the  compres¬ 
sion  filter  appear  to  be  at  least  23  db  lower  in  amplitude  than  the  com¬ 
pressed  pulse. 

C.  CONCLUSIONS  FROM  EXPERIMENTAL  RESULTS 

The  test  filter  performance  provided  a  substantial  improvement  over 
that  of  other  techniques  previously  used  at  Stanford  to  obtain  pulse  com¬ 
pression  at  microwave  frequencies.  The  most  serious  drawback  of  this 
type  of  microwave  filter  is  the  relatively  high  insertion  loss.  The  use 
of  directional  couplers  as  taps  would  alleviate  this  objection  while 
simultaneously  providing  an  easily  adjustable  weighting.  It  appears 
that  a  filter  with  a  considerably  larger  time-bandwidth  product  could  be 
constructed  using  this  technique  with  a  longer  delay  line  and  more  taps, 
the  limiting  factor  probably  being  the  accuracy  with  which  the  excitation 
function  could  be  generated. 
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IV.  CONCLUSIONS  AND  REOOIMENDATIONS 


Hie  computational  results  described  indicate  that  the  tapped-delay- 
line  compression  filter  can  perform  competitively  with  many  other  systems. 
Some  systems  making  use  of  tapped  delay  lines  have  used  bandpass  filters 
on  the  individual  taps  in  order  to  suppress  the  grating  lobes  that  are 
produced  as  a  result  of  the  low  .number  of  taps.  It  is  possible  that  this 
apparent  preference  for  the  tuned -tap  compression  filter  is  due  to  the 
difficulty  of  predicting  the  response  of  the  untuned- tap  filter  by  ana¬ 
lytic  methods.  The  advent  of  the  high-speed  computer  has  eliminated  this 
objection  to  the  use  of  the  untuned- tap  filter  and  should  allow  this 
filter  to  be  seriously  considered  for  many  pulse-compression  applications. 

The  experimental  filter  has  shown  that  even  a  very  simple  construc¬ 
tion  technique  can  provide  a  compression  filter  whose  performance  rivals 
much  more  complicated  systems  [Ref.  4].  No  attempt  was  made  with  the 
experimental  filter  to  realize  the  potential  performance  of  the  tapped 
delay  line  due  to  the  limitations  imposed  upon  the  compression  system  by 
other  components. 

The  recent  development  of  much  improved  UNO's  makes  it  possible  to 
generate  broadband  (4  Gc)  frequency  sweeps  with  a  maximum  frequency 
deviation  from  linearity  of  only  10.1  percent  of  the  total  bandwidth 
swept.  Measurements  on  one  tube*  indicate  that  this  linearity  can  be 
maintained  at  scan  rates  greater  than  20  Mc/nsec.  The  fine-grain  power 
variations  of  this  tube  are  also  much  improved  over  that  of  the  system 
used  for  testing  the  experimental  filter.  Reduction  of  these  power 
variations  would  have  the  effect  of  reducing  the  spurious  responses 
present  in  the  compression-filter  output  due  to  direct  detection  of  the 
local-oscillator  rf  envelope. 

These  improvements  make  it  possible  to  realize  more  of  the  potential 
of  the  compression  filter  in  an  actual  system.  Thus  it  appears  that  the 
construction  of  a  filter  with  the  weightings  derived  by  the  computer 
program  would  justify  the  effort  if  the  new  filter  were  used  in  a  system 

« 

A  Varian  VA-175G  backward  wave  oscillator  manufactured  by  Varian 
Associates,  Palo  Alto,  California. 
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employing  presently  available  BWO's.  Such  an  improved  system  might  still 
be  limited  by  internal  reflections  in  the  TWT  amplifiers,  although  not  to 
the  extent  that  the  test  filter  described  here  was  limited.  Because  of 
the  hirgh  losses  in  the  experimental  filter,  the  TWT  amplifier  had  to  be 
driven  to  saturation  in  order  to  obtain  a  sufficient  signal-to-noise 
ratio  at  the  sampling  oscilloscope  for  observing  the  sidelobes.  The 
spurious  signals  generated  by  reflection  within  the  TWT  appear  about  20 
db  lower  in  amplitude  than  the  primary  signal  when  the  tube  is  near 
saturation,  but  drop  an  additional  10  db  lower  if  the  power  level  is 
kept  at  least  10  db  below  saturation.  Several  manufacturers  have  recently 
placed  tubes  on  the  market  which  are  capable  of  delivering  10  w  over  an 
octave  band  while  maintaining  noise  figures  of  15  db  or  better.  One  of 
these  tubes  should  nearly  eliminate  the  TWT  as  a  source  of  sidelobes. 

The  tapped  delay  line,  although  treated  in  this  report  as  a  micro- 
wave  device,  is  certainly  not  restricted  to  microwave  frequencies.  A 
100-tap-delay-line  compression  filter  centered  at  500  kc  has  been  used 
in  a  scanning  spectrum  analyzer  system  at  Stanford  [Ref.  5].  Corning 
Glass  Works  [Ref.  6]  describes  a  glass  ultrasonic  delay  line  with  as  many 
as  1000  taps  and  a  potential  bandwidth  of  about  20  Me.  The  large  numbers 
of  taps  involved  with  these  filters  make  the  contribution  of  any  individ¬ 
ual  tap  a  very  small  percentage  of  the  total  response.  The  computer  then 
becomes  the  only  practical  method  of  finding  the  tap  weighting  which  pro¬ 
vides  the  best  response  for  a  given  application. 
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APPENDIX  A.  COMPUTER -PROGRAM  FLOW  CHART 


1 .  Def ini tions 
A 

A1 

A2 

AMPC 

AMPLNG 

AMPRNG 

AP 

ARNGID 

AT 

ATICK() 

a(tmtau) 

axplot(arng) 

B1 
B2 
BP( ) 
bt(  ) 
btick() 


of  Terms  Used  in  Flow  Chart 

if  -10  >  TWTAU  >60,  A  =  0.0  else  A  =  1,  or 
A  =  1  -  [0.0081632653  (TOTAU  -  25)2] 

minimum  value  of  amplitude  range 

maximum  value  of  amplitude  range 

amount  of  amplitude  equal  to  1  in.  on  amplitude  axis 

length  of  amplitude  axis  in  inches  (ilO.O) 

reciprocal  of  amplitude  range  (scale  factor  for  total 
plot) 

parameter  of  PKCALC  used  to  obtain  valid  set  of  t's 
or  B's  depending  upon  which  is  being  perturbed 

scale  factor  of  amplitude  for  plot  which  skips  the 
main  lobe 

parameter  of  PKCALC  used  to  obtain  valid  set  of  T's 
or  B's  depending  upon  which  is  being  perturbed  at  a 
given  time 

values  of  10,  5,  2.5,  1.25  or  1  (in  percent)  of  PKMAX1 
used  for  tick  marks  on  left  side  of  plot 

calculate  multiplier  of  cos  function  ESUM 

plots  time  axis  and  box 

current  value  of  BP() 

current  value  of  BT()  if  perturbing  B,  else  holds  B1 
matrix  of  valid  set  of  B  values 
matrix  of  tested  B  values 

same  as  ATICK()  except  for  right  side  of  plot 
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CHANGE 

CHNG 

E0(t) 

EREV 

ESUM 

esum(time) 

ET1 
ET2 
ET3 
ET4 
IB 

ICNT 
I  FLAG 

IM1 
IPERM 
ITAU 
I  TEMP 
IT12 
ITM1 

11M1M1 

i 

I1M2 

I 
! 
i 


decimal  representation  of  percent  change  peak-to- 
sidelobe  ratio  for  cutoff 

change  in  peak- to-sidelobe  ratio  magnitude 

value  of  r(t)  with  second  index  used  to  indicate 
permanent  or  temporary  set--thereby  used  as  a  switch 
to  indicate  valid  function 

calculate  r(t)  for  perturbed  T  or  B 

r(t) 

calculate  r(t)  at  T  =  TIME  for  given  TAUP()  and  BP() 
(initial  summation) 

r(t)  of  first  point  of  main  lobe,  valid 
r(t)  of  last  point  of  main  lobe,  valid 
scaled  value  of  ET1 
scaled  value  of  ET2 

flag  for  calculating  or  reading  in  B's 

counts  the  number  of  perturbations  of  T's  and  B's 

flag  to  indicate  which  set  is  being  perturbed  (o  for 
T,  1  for  B) 

I  minus  1 

index  to  indicate  valid  set  of  given  matrix  quantity 

flag  for  calculating  or  reading  in  t's 

index  to  indicate  tested  set  of  given  matrix  quantity 

number  of  time  points  within  area  of  interest 

matrix  value  of  first  point  to  begin  expecting  main 
lobe 

last  matrix  point  before  main  lobe  section 
matrix  value  of  last  point  to  expect  main  lobe 
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first  matrix  point  after  main  lobe  section 


LASTT 

number  of  tick  marks  on  time  axis 

MAXCNT 

number  of  runs  through  T's  and  B's  allowed 

N1 

first  tap  number 

N1NC 

increment  of  n 

N1P1 

second  tap 

N2 

last  tap  number  ??  101 

NOPKS=NPK 

the  quantity  on  the  left  is  the  temporary  equivalent 

of  the  one  on  the  right  within  PKCALC  when  a  new  set 

of  T's  or  B's  is  being  tried.  Changing  a  set  in 

value  replaces  the  set  on  the  right  by  that  on  the 

left 

NPK 

number  of  peaks 

NPK1 

number  of  peaks  in  area  before  main  lobe  is  encountered 

NPK2 

number  of  peaks  from  start  through  end  of  main  lobe 

area 

N2M1 

next  to  last  tap 

pkcalc(ap.at) 

calculate  new  possible  set  of  peaks 

PKMAG ( , ) 

matrix  to  hold  peak  values--last  index  to  switch  valid 

groups 

PKMAX1 

magnitude  of  maximum  peak- -eventually  main  lobe 

PKMAX2 

magnitude  of  second  largest  peak- -eventually  maximum 

sidelobe 

pkrat(  ) 

matrix  to  hold  ratio  of  main  lobe  to  peaks 

PKTIM1 

time  of  PKMAX1 

PKTIM2 

time  of  PKMAX2 

PL0T1 

plots  initial  and  final  plots 
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PIjOT2 

PRES 

PREV 

PROD 

RATIO 

RATIOl 

SAVE 

SPKMAX 

S  TPMAX =SPKMAX 
T() 

TAU 

TAU1 

TAU2 

tau(fln) 

TAUINC 
TAUMAX 
TAUP( ) 

TAUT() 

TICKS 

TIME 

TINC 

TL 


plots  intermediate  data 

indicator  by  sign  of  slope  between  last  two  points 

indicator  by  sign  of  slope  between  the  previous  last 
two  points 

indicator  of  whether  a  perturbation  will  increase  or 
decrease  the  value  of  T 

ratio  at  beginning  of  a  set  of  perturbations 

current  ratio  for  best  case  (peak  main  lobe/maximum 
sidelobe) 

holds  the  last  magnitude  of  the  last  peak 

actual  value  of  PKMAX1 

same  as  NOPKS=NPK 

matrix  of  time  value 

100  /l. 01  -  0.0l(N)  -  50.0 

current  value  of  TAUP() 

current  value  of  TAUT()  if  perturbing  T,  else  holds 
TAU1 

calculate  TAUP()  if  not  read  in 

perturbation  increment  of  T 

maximum  allowable  value  of  any  T 

matrix  of  valid  set  of  T  values 

matrix  of  tested  set  of  T  values 

calculates  the  values  of  the  tick  marks  on  the 
amplitude  axis 

T(),  for  calculation  of  r(t) 

time  increment  for  all  time' variables  (sec) 

reciprocal  of  length  of  time  axis  (used  in  plotting) 
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TL1 


v 


TLAG 
TLNGTH 
TWA INI 
TMAIN2 
TWIN  I  ,TMAXI 

TTHTAU 

TMTAU1 

TMTAU2 

TPMAX1=PKMAX1 
TPMAX2 =PKMAX2 
TPTIM1=PKTIM1 
TPTIM2=PKTIM2 
TR 

TRATI0=RATI01 

TRNGE 

TZERO 


length  of  time  axis  plus  0.05  in.  (used  in  plot  ticks) 

plotting  constant  for  time  tick  marks 

length  of  plot-time  axis  (inches) 

minimum  time  to  begin  expecting  main  lobe 

time  to  stop  expecting  main  lobe 

beginning  and  end  of  the  time  interval  within  which 
the  peak- to-sidelobe  ratio  is  optimized 

T  -  r 

time  minus  TAU1 
time  minus  TAU2 
same  as  N0PKS=NPK 
same  as  N0PKS=NPK 
same  as  N0PKS=NPK 
same  as  NOPKS  xNPK 
total  range  of  time 
same  as  N0PKS=NPK 

amount  of  time  equivalent  to  1  in.  on  time  axis 

number  of  time- axis  inches  equal  to  the  minimum  time 
used  to  shift  time  axis  to  zero 


2.  Flow  Chart  Diagram 

A  detailed  flow  chart,  which  describes  in  nonsyntactical  language 
the  program  listed  in  Appendix  B,  is  given  on  pages  40  through  58. 
Descriptive  notations,  as  well  as  statement  numbers,  are  included  within 
the  flow  chart  where  appropriate. 
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DETAILED  FLOW  CHART  OF  COMPUTER  PROGRAMS  AND  SUBPROGRAMS 


Set  arrays  T(250l\  TAUP(lOl',  TAUT(lOl\ 
BP(lOl),  BT(lOl',  E012501 ,2  V,  PKMAG ( 2000 . 2 ' . 
PKTIM(  2000 , 2) ,  PKRAT(  2000 ' ,  ATICK(20'I.  BTICK(,20' 


Read  In  parameters  /  Read  in 

concerning  time  and  (  TMINI ,  TMAXI,  TINC,  TMAIN1, 
plot.  \  TMAIN2,  TLN’GTH,  AMPLNG 


Calculate  variables  for 
plot . 


TR  =  TMAXI  -  TMINI 

1 

TL  =  1.0/TLNGTH 

* 

TL1  =  TLNGTH  +0.05 

\ 

TRNGE  =  1.0/TR 

t 

TZERO  =  |  TMINI  j  \  TRNGE 

t 

AMPLNG 

A1  =  5.0  -  2 

1 

A2  =  A1  +  AMPLNG 

* 

AMPC  =  0.05  X  AMPLNG 

1 

mo  »  -tlag 

%.  TTAG  :  0 y 

TLAG  =  10.0  -  HAG 
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m 


LASTT  =  0.1TR  +  1.05 


ITI2  =  TR/TINC  +1.5 


Calculate  number 
of  time  intervals 


TMAIN1  -  TMINI  .  . 
1  =  - TINC -  +  15 


Time  interval  at 
start  of  main  lobe. 


TMAIN2  -  TOIN1  .  _ 

2  _  - — — -  +  1.5 

TINC 


rval  at 
ain  lobe. 


Tj  =  THIN1 


Generate  the  times 
for  each  time 
interval . 


I  =  2 


Tj  =  Ti  ♦  TINC 


I  s  I  »  1 


I  TNI Ml  =  I TNI  -  1 

1 

1TM2P1  «  IHI2  +  1 

t - 

ICNT  «  0 

Last  time  interval 
before  main  lobe. 


First  time  interval 
after  main  lobe. 

Initialize  the 
number  of  times  t 
and  B  perturbed . 
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IB  :  0 
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Sense  light  4  is  a 
switch  indicating 
whether  working  on 
area  before  main 
lobe  or  area  after 
main  lobe. 


i 


53 

PRES  = 

E02,l 

'  E°l,l 

SENSE 

LIGHT 

4  ON 

Slope  between  first  two 
points.  It  will  be  used 
to  determine  if  a  peak 
has  been  reached. 


Initialize  II  to  3 
because  EO  has  already 
been  found  for  T1  and 

T2 

Set  12  to  last  T  incre¬ 
ment  before  main  lobe. 


This  loop  calculates 
EO  before  and  after 
main  lobe. 

Set  PREV  to  the 
slope  just  found. 


Calculate  EO  for 
next  T  interval. 


Find  slope  of  this 
time  interval. 
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If  the  slope  is  unchanged, 
there  is  no  peak. 


If  slope  changes,  the 
last  EO  was  a  peak. 


Store  magnitude 
of  peak. 


Store  time  of  peak. 


Use  SAVE  to  test 
for  largest  peak. 


If  SAVE  >  | value |  of 
the  largest  peak,  make 
it  the  largest  peak. 
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Have  finished  EO  with 
the  initial  conditions 


I 


I 


Must  do  the  main  lobe 
and  area  after  the  main 
lobe. 


Store  the  value  of  the 
EO  encountered  just 
before  the  main  lobe. 


Store  the  magnitude  of 
this  point. 


Store  the  time  of  this 
point. 


Store  the  number  of 
peaks  before  the  main 
lobe.  This  will  be  used 
to  plot  the  peaks  with¬ 
out  the  main  lobe. 


This  loop  calculates  EO 
of  the  main  lobe  and 
finds  the  value  and  time 
of  the  main  lobe. 
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m 


No  change  in  slope 
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Write  out  the  initial 
calculations . 


Main  and  sidelobe 
peak  information. 


The  values  of  B  used. 


The  values  of  t  used . 


Calculate  the  tick  marks 
for  the  main  plot . 


Plots  EO  vs  time.  Also 
plots  the  peak  values 
without  the  main  lobe. 
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RATIO  holds  the  ratio  of 
the  last  set  of  B's  and  T's 
before  this  perturbation. 


Increment  the  number  of 
times  B  and  T  perturbed. 


Switch  indicating  B's  have 
been  perturbed  but  not  t's, 


r 

l 
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This  section  perturbs 
the  T's,  tests  them 
and  plots  intermediate 
plot.  If  within 
change  requested, 
plots  final  plot  and 
stops . 


© 


RATIO  =  RATI01 


I  FLAG  =  0 

J1  = 

N1 

Save  new  ratio. 


Now  perturbing  the  T's. 

Do  the  first  tap  outside 
the  loop  since  may  have 
to  check  against  thu  min 
and  max  T. 


TAUT,,  =  TAUP,,  +  TAUINC 
J 1 


90 


100 


<sTAUTt1  :  TAUMIN 

1 

95 

TAUTN1  not  >  TAUMAX 

105 

TAUT,,  =  TAUMIN 

or  <  TAUMIN 

TAUT,,  =  TAUMAX 

J1 

Calculate  a  new  set  of 
peaks  using  new  TAUTj^. 

If  better,  PKCALC  will 
make  them  a  valid  set. 


r  110  ' 

PKCALC  (TAUP,,  ,  TAUT..) 

L  J1  I 


The  original  TAU  was 
better.  Try  subtracting 
TAUINC . 


The  TAUTN1  was  better 
than  TAUPNj^  ,  changes  have 
been  made;  go  to  the  rest 
of  the  TAUT. 


SENSE  LIGHT  3  OFF 
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22 


New  TAUTjj^  does  not 
need  to  be  checked . 


.. ,  1 

140 

Check  TAUTj^  so 

- 1 - 

130 

TAUTt1  =  TAUMIN 

J  1 

that  not  >  TAUMAX 

TAUTt1  =  TAUMAX 

or  <  TAUMIN 

145 

PKCALC(tAUPt, ,TAUTt1 ) 
J 1 _ J  J. 


0- 


1 

J  = 

50 

N1P1 

J1 

=  J 

-© 


Generate  new  peaks. 
If  better,  PKCALC 
will  make  changes. 


TAUT, 

=  TAUP,  +  TAUINC 

J 

J 

^  PKCALC (TAUPj, TAUT j)  ^ 


This  loop  perturbs  the 
T's  from  the  second  to 
the  next  to  last. 

These  are  not  checked 
against  TAUMIN  or 
TAUMAX,  but  the  method 
is  the  same. 


£ENSE  LIGHT  J  ON 


155 

TAUT.  =  TAUP.  -  TAUINC 

J  J 

f  " 

^  pkcalc(tau 

PJ<  tautj)  ) 

u 

1 

J  -  J 

♦  NINC 

J  :  N2MI 


J1  =  N2 


This  section  perturbs  the 
last  T.  Checks  against 
TAUMAX  and  TAUMIN  if 
necessary . 


TAUT  =  TAUPt1  +  TAUINC 

J  1  «J1 


PROD  :  0 


- - -  ^ 

TAUINC  :  0 


5^UTJ1  :  TAUMIN 


NO 


26 
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Write  out  the  last 
change  in  ratio  and 
which  set  was  last 
to  be  perturbed, 


300 

WRITE  OUT  (PRINTER) 

"IN  THE  "  ICNT ,  "30TH  RUN,  THE  CHANGE 
IN  RATIO  =  "CHNG , "  WITH  THE  LAST 
PERTURBED  SET  BEING  "  I FLAG 


WRITE  OUT  (PRINTER) 

"MAIN  LOBE  MAG.  =  "PKMAX1 
"AT  T  ="  PKTIM1 
SKIP  A  LINE 

"MAX.  SI DELOBE  MAG.  ="  PKMAX2 
"AT  T  ="  PKTIM2 
"GIVING  A  RATIO  »  RATIOl 


Print  out  main  and 
sidelobe  information 
and  ratio. 


WRITE  OUT  (PRINTER) 
"THE  VALUES  OF  B(n)" 
Bn  N  =  Nl,  N2,  NINC 


Print  out  the  values 
of  B. 


WRITE  OUT  (PRINTER) 
"THE  VALUES  OF  TAU(N)" 
TAUPj  I  =  Nl,  N2,  NINC 


Print  the  values 
of  T. 


Plots  the  last  set 
of  EO's. 
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APPENDIX  B. 


FORTRAN  PROGRAM  LISTINGS 


1 .  Input-Output  Data 

The  input  data  for  the  FORTRAN  program  are  tabulated 


Card  Location 


Number  (Columns) 

Variable 

Definition  of  Variable 

1  1-8 

TTAINI 

Time  to  begin  testing 
peak-to-sidelobe  ratio 

9-16 

TMAXI 

Time  to  stop  testing 
peak-to-sidelobe  ratio 

17-24 

TINC 

Time  increment  for  all 

time  variables 

25-32 

TMINC 

Initial  value  of  time 
for  plot  purposes 

33-40 

TMAXC 

Final  value  of  time 
for  plot  purposes 

41-48 

TMINP 

Initial  value  of  plot 
time  axis 

2  1-8 

TMAXP 

Final  value  of  plot 
time  axis 

9-16 

TLNGTH 

Length  of  plot  time 
axis  (inches) 

17-24 

TOAIN1 

Beginning  of  time 
interval  within  which 
main  lobe  is  deter¬ 
mined 

25-32 

TMAIN2 

End  of  time  interval 
begun  by  TMAIN1 

33-40 

AMPLING 

Height  of  amplitude 
axis  in  inches 
(510.00) 

below. 


Format 


+  XXXX . XX 


ff 


It 


II 


It 


If 


It 


tl 


II 


II 


It 
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Card  Location 


Number  (Columns) 

Variable 

Definition  of  Variable 

Format 

3  1-5 

N1 

Value  of  first  n 
(integer) 

xxxxx 

6  10 

N2 

Value  of  last  n 
(integer) 

IT 

11-15 

NINC 

Spacing  between  n's 
(integer) 

It 

16-20 

ITAU 

Boolean  variable  (0  if 
initial  T's  are  cal¬ 
culated,  1  if  they  are 
read  as  data) 

tt 

21-25 

IB 

Boolean  variable  (0  if 
initial  Bn's  are  0.5, 

1  if  they  are  read  in) 

ft 

26-30 

MAXCNT 

Maximum  number  of 
times  to  run  through 

V  Bn  set 

tt 

4  1-14 

TAUINC 

Perturbation  increment 
for  T's 

±0.xxxxxxE±xx 

15-28 

BINC 

Perturbation  increment 
for  B's 

tt 

29-42 

TAUMIN 

Minimum  allowable 

value  for  T 

tt 

43-56 

TAUMAX 

Maximum  allowable 
value  for  T 

tt 

57-70 

CHANGE 

Change  in  peak-to- 
sidelobe  ratio  below 
which  program  will 
terminate 

tt 

5  1-14 

TAUP(l) 

First  T  (if  ITAU  = l) 

tt 

(The  T's  are  arranged  in  increasing  order  of  subscripts,  five 
to  a  card,  on  as  many  cards  as  necessary  to  contain  them.  They 
K  are  contained  in  columns  1-14,  15-28,  29-42,  43-56,  and  57-70.) 


K+l  1-6  BP(l)  First  B  (if  IB  =  l)  +xxx.xxx 

(The  B's  are  arranged  in  increasing  order  of  subscripts,  ten 
to  a  card,  in  columns  1-6,  7-13,  14-20,  21-27,  28-34,  35-41, 

Z  42-48,  49-55,  56-62,  and  63-69.) 
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An  outline  of  the  output  data  for  the  FORTRAN  program  is  given 


below. 


a.  Initial  Case 

(1)  Printed  Output 

(a)  Magnitude  and  time  location  of  main  lobe. 

(b)  Magnitude  and  time  location  of  largest  sidelobe. 

(c)  Initial  values  of  T's. 

(d)  Initial  values  of  B's. 

(2)  Plotted  Output 

(a)  The  function  r(t)  is  plotted  from  H4INP  to  TMAXP 
in  intervals  of  TINC.  The  plot  is  contained 
within  a  box  of  length  TLNGTH  and  height  AMPLNG. 

The  height  of  the  main  peak  is  scaled  to  4.5  in. 

(b)  The  function  r(t)  is  plotted  as  above  except  the 
amplitude  scale  is  increased  X10  and  the  section 
from  1MAIN1  to  TMAIN2  is  omitted. 

b.  Intermediate  and  Final  Cases 

(1)  Printed  Output 

Same  as  initial  except  that  the  increase  in  peak-to-sidelobe 

ratio  over  the  previous  case  is  printed,  and  only  the  param¬ 
eters  (t's  or  B's)  perturbed  on  the  last  case  are 

printed . 

(2)  Plotted  Output 

Same  as  initial  case. 


2.  Program  I 

This  program  calculates  the  function 

m 

r(t)  =  )  A(t-T  )  B  cos  2jt  2(t-T  )  -  O.Ol(t-T  )2  +  <> 

^  n  n  L  n  n 

n=l 

>  for  various  values  of  $.  All  values  of  B  and  T  must  be  read  in 

n  n 

as  data,  and  the  different  values  of  4  must  be  inserted  in  the  body  of 
the  program  as  PHl().  No  flow  chart  is  provided  for  this  program  due 
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to  its  similarity  to  Program  II.  No  perturbations  are  done--the  function 
is  calculated  once  and  plotted.  Subprograms  required  are  ESUMII, 

A(t-TAU) ,  PL0T1 ,  AXPLOT,  TICKS,  as  well  as  the  three  Stanford  plot 
routines  XYPLOT,  XYLABL,  and  XYFRC . 


01HENSI0W  Tc250l)»TAUP(101),TAUTC101),aP<101>,BT(101)»£0<2501»2>» 
1PKMAG{2000,2),P<TIM(2000»2).PKRATC2000), AT  I CKC 20 ) * BT I CK C 20 ) , 
2PHK21) 

COMMON  TWAXI*T1NC»TR»TLNGTh»TRNGE»TLAG»LASTT,IT12»ITM1*ITM2, 
1ITM1M1,ITM2P1,N1*N2.NINC»NP«,NPw1,NP<2,P<MAX1,Pi<MAX2,P«TIM1, 
2PKTIM2>1oERM«ITEnP,U1»IFLAG»AMPRNG»A1»A2»ARNG10»RATI01,5PWMAX» 
3TZE«0#T»TAJP»TAJT.B»»BT,EO#P'<MAG»PKT|M»P<BAT*ATICM#BT1CK#TL#TL1» 
4PHI,PHII 

C  RE A0  In  basic  time  PARAMETERS.  SET  up  PLOTS  AND  INDEX  PARAMETERS 

READ  INPUT  TAPE  5. 500,  TM I NI , TM AX  I , T I NC . TMAI Ni  ,  TMA I N2, TLNGTH, AMPLNG 

tr-tmaxi-tmini 

TL-1.0/TLN3TH 
TL1»TlNGTH»0.05 
TRNGE-1 ,0/TR 
T2ERO«A0SF(TMINI }*TRNGE 
A1»5.0-0.5*AMPLNG 
A2-A1*AMPlNG 
AMPC»0,05*AMPlNG 
TLAG-M0DFCTMINI.10.05 
If(TLAG)l»10»5 
1  TLAG— TLA3 
GO  TO  10 

5  tlag-io.o-tlag 
10  LASTT-0.1-TRM.05 
I TI 2-TR/T I NC* 1 • 5 
ITM1-CTMAIN1-TMIND/TINC41.5 
ITM2-CTmAIN2-TMINI)/TINC*1.5 
TCD-TMINI 
00  15  1-2. 1  TI 2 
15  TCI)-TU-1)*TINC 

ITM1M1-ITM1-1 

ITM2P1-ITM2*1 

ICNT-0 

C  READ  IN  REMAINING  PARAMETERS,  SET  UP  TAU  AND  8 
REAO  INPUT  TAPE  5, 510, NI , N2, NI NC 
N1P1-N1+NINC 
N2M1-N2-1 

20  REAO  INPUT  TAPE  5,520, C TAUPC I >» I-Ni, N2, NINC > 

AO  REAO  INPUT  TAPE  5,530,  ( BP( I), I >Nl ,N2,NINC ) 

WRITE  OUTPUT  TAPE  6, 620,  C BP{ I > , I -Ni , N2, NINC 5 

WRITE  OUTPUT  TAPE  6, 6 1 0,  < T AUP( I ) » I -N1 » ' . . NI NC 5 

NPH I  ■  10 

PH  I ( 1  )  -  0,10480 

PHI <  2  )  ■  0,15011 

PHI(S)  -  0,01536 

PHI C 4 )  -  0,02011 

PHI ( 5  5  -  0,81647 

PH X  C  A  )  ■  0,91646 

PHICTJ  ■  0,69179 

PHI(S)  -  0,14194 

PHIC9)  -  0,62590 

PHI ( 10 )  ■  0,36207 

c  calculate  initial  sums,  generate  initial  plots  and  side  lobe  structure 

00  315  K-1,NPHI 


PHI  I  •  PHI(K) 

47  NPK*0 

eou»i)*esjm<tci>> 

E0(2»1)*ESJH(T(2)) 

IF(A8SF(EO(2#n)-A0SF(EOCl»l))>52»51»51 

51  PKMAX2*A9Se,(E0{2»ll) 

PKTIM2«T(2) 

GO  TO  53 

52  PKMAX2*A3SrtE0<l#l)> 

PKTIM2«T(1) 

53  PRES«E0(2#n-E0(l*l) 

SENSE  LIGHT  4 

11*3 

1 2* I TM 1  Ml 

54  00  75  1*11.12 
I  Ml ■ I • 1 
PREV-PRES 
E0(I»1)»ESJM(T(I)) 

PRES«EO(I#n-EO(IMl#l) 

IF(PR£V*PRES)55.55.75 

55  NP<*NPK*1 
PKMAG(NP<»1)*E0CIM1.1) 

pktim(np<*d«t(imi> 

SAVE*A3SFC3<MAG<NP<»1>> 

IF(5A9E*pKHAX2)75#70»70 
70  PKMAX2*SAV£ 

PKTIM2-TCH1) 

75  CONTINUE 

IF{ SENSE  LIGHT  4)77.02 

77  SP<MAX*E0<ITM1M1.1) 

PKMAXl«A3Sf(SPKMAX) 

PKTIM1«T(ITM1M1) 

NPXl-NPK 

00  80  I«IT41#ITH2 

IM1-I-1 

PREV*PRES 

E0( I. 1 )*ESJH(T( I ) ) 

PRES»E0(I»1)-E0CIM1,1) 

IF(PREV*pRE5)70*78.3O 

78  NPK*NP<*1 
PKMAG<NP<»l)»EO(IMl.t> 

PKTIM(NP<#l)«TCIMl) 

SAVE»ABSF(E0<IM1*1)> 

IF(SAVE“P|<HAX1)80»79»79 

79  PKMAX1-SAVE 
PKTIMl»T(Hl) 

SP<MAX»EO(IMi#l) 

00  CONTINUE 
I1»ITM2P1 
1 2*1 T I  2 
NPK2*NPKM 
GO  TO  54 
82  IPERM-1 
ITEMP*2 

PATI01-PXMAX1/PXNAX2 

AMPRNG-AMPC/PXMAX1 

IF( RAT  101-10.0)1000. 1001* 1001 

1000  ARNG10-AMPRNG 
GO  TO  1002 

1001  ARNG10*10,0*AMPRNG 

1002  HRITE  OUTPJT  TAPE  6.605.PKMAX1.PKTIM1.PKMAX2.PKTIM2.RATI01 
CALL  ticxs 

CALL  PLOTl 
315  CONTINUE 
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500  F0RMATC7F8.2) 

510  F0RMAT<315> 

520  F0RMATC5E14.6) 

530  FORMAT ( 1 0F6 , 3 ) 

605  FORMAT ( 18H0MAJN  LOBE  M AG .  ■  # E 15. 8, 2X» 7HAT  T  ■  ,FB. 3/23H0MAX.  SIDE 
1  LOBE  HAS.  «  ,E15.8,2X,7HAT  T  ■  , FB , 3, 5X, 17HGI VING  A  RATIO  «  »E15, 
28) 

610  F0RMATC21H0THE  VALUES  OF  TAU(N)/C10E13,6) ) 

620  F0RMAT(19H0TH£  VALUES  OF  B< N)/ C 20F6 . 3 )  ) 

625  FORMAT (14) 

630  FORMAT(lOF8.5) 

CALL  EXIT 
END 


3.  Program  II 

This  is  the  general  peak- to-sidelobe  ratio  optimizing  program.  The 
function  calculated  is 


m 

2(t-Tn)  -  0  Ol(t-Tn)2] 

n=l 


'(t)  =  ^  )  B^  cos  2ji 


The  initial  values  of  B  must  be  read  in  as  data,  and  the  initial 

n 


values  of  T  are  calculated  by  Subprogram  TAUN.  The  region  of  the 
n 

main  peak  is  TMAIN1  <  t  <  TMAIN2,  and  the  sidelobe  is  considered  the 
largest  peak  outside  this  region.  This  program  will  terminate  when  a 
series  of  perturbations  results  in  an  improvement  of  the  peak-to-sidelobe 
ratio  less  than  CHANGE.  Subprograms  required  are  ESUMI ,  A(t-TAU),  PIX3T1 , 
PLOT2,  AXPLOT,  TICKS,  XYPLOT,  XYLABL ,  and  XYFRC. 


DIMENSION  T  I  2*iC  1 )  ,TAUP<  10 U  ,TAUT(  101*  ,HP(  lol)  ,I>T(  1C1 1  •  liO (  2501,2)  , 
1PKFAU  (  2000. 2  ) ,  PK  T I M  (  2000 , 2  ) ,  PKP.  AT  ( 2000 ) ,  AT  ICK  l  20 )  , bT ICK ( 20 ) 

COMMON  TMAXl  .TWC.TR,  r  LNUTh ,  TKNGfc ,  T  L  AG ,  L  AST  r ,  I  T  1 2 , 1  ftol , ITM2 , 

1 1 TH1M1 » I  rM2Pl,M,N?,NlNC»f|PK,NPKl,NPK2,PKMAXl,PKMAX2,PKT  IM1, 
2PKTIM2, 1  PERM,  1 IEMP,  J1 ,  IFIAG,AMPRNG,  Al,A2,  AR.'UilC',KAT  IUl.SPKtiAX, 

3  TZEKO , 1 »  T  AUP ,  TAtll  ,  HP , « F , 10.PKMAG, PKTIM, PKRAT ,  A  I  ICK ,  iH  ICK,  1 L  ,  1  L 1 
READ  IN  BASIC  1  1  ME  PARAMETERS,  SET  UP  PLOTS  An'.)  INDEX  PAkAMETEk., 

READ  INPUT  TAPE  5 , 500 ,  TrtI  N I ,  TMAX I ,  TINC ,  IMA  I  4 1 , 1  MA I N2 ,  TLNG  TH ,  AMPL.NG 

TR=TMAXl-TMiNI 

TL=1.0/TLNGTH 

TLl*TL.NGTH«-O.Or, 

TRNGE-l.O/TK 
T4!tl<0=ABSF(  TMINl  )  *TKNG£ 

AI=5.0-0.5*AMPLNG 
A2=AUA“PLNG 
AMPC-0. )5*AMPLN0 
TLAG  =  MUDF( TMINl , 10. 0) 

I F ( FLAG) l, lu,5 
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1  TLA'j=— TLAG 
GO  TO  10 

5  T L  AG  =  LC . 0-1  LAG 
1C  LASTT  =  0.  l*Tm-l .05 
I T 1 2=  TK / T I NC ♦ 1  •  5 

I  T  M 1  =  (  I  HA I N 1-  IN  N I )  /  r  I NC  +  1 . 5 
I TM2= (  TMAIN2-1>’INI)/T  1NC+1 .5 
T (  1)=TM1NI 

DC  1K  I=2,ITI2 
15  r<  1  )=T(  I-1J+TIKC 
I IKlKl=i TM1-1 
ITP2P1=I  IM2+1 
I C  M  =  0 

C  KLAD  li\  REMAINING  PAR AML  f  t K 3 ,  SET  UP  TAU  AND  t 

I’C  All  INPUT  TAPL  5 , 5 10  ,  N 1 ,  N2 ,  N  INC .  I  TAU,  1 H ,  MAXLN 1 ,  T  AD  I NC  ,  B I NC  ,  T  AUM  IN 
1 , T  AUM AX , CHANGE 
N1 P 1 =N1+NI NC 
N2M1=N2-1 
1 F  I  1TAU120, 25,20 

20  READ  INPUT  IAPl  5 , 5 20 , t T AUP (I) , 1 =N 1 ,N2 , N1 NC ) 

GC  rc  35 

25  DC  30  I =N1 , N2 , M NC 
FL  1  =  I 

TAUIM  I )  =  TAU(FL I  I 
30  CONTINUE 

35  PKCO=TAUINC*( T AUP (N2 1 -T AUP 1 Nl ) ) 

I F  C IB  >40,45,40 

40  READ  INPUT  TAPE  5 , 5 30 , ( HP (I) , I =  N1 , N2, N I NC ) 

GO  TO  47 

45  DO  4ft  I  =N  1 ,  !'<2 , NINC 

46  BP ( I)=0.5 

C  CALCULATE  INITIAL  SUMS,  GENEKATL  INITIAL  PLOTS  AND  SIDE  LOBE  STRUCTURE 

47  NPK  =0 

FO (  1 , 1 )  =E$UM( T (  1 1) 

EO  (  2, 1 )  =ESliM  (  T  (  2 )  ) 

I F ( ArtSF ( EO ( 2 , 1 ) l-AHSF ( E0( 1,1)1)52,51,51 

51  PKMAX2  =  AHSF(EO(2,l)  ) 

PKTIM2=T(2) 

GO  TO  53 

52  PKMAX2=ABSF(E0( 1,  l)  ) 

PK  T I M2=  T ( 1 ) 

53  PRES=EC(2, 1 )-E0( 1,1 ) 

SENSF  LIGHT  4 

I I  =  3 

I2=ITM1M1 

54  00  75  1=11,12 
IH  1  =  1-1 
PREV=PXES 

EO ( 1 , 1 ) =ESUM( T ( I )  ) 

PRcS=EO(I,l)-tC«IMl,l) 

IF(PREV*PRES)55,55,75 

55  NPK=NPK+ 1 
PKNAG(NPK,l)=EO(IMl»l) 

PK  T I M<  NPK , 1 )  =  T ( IMl) 

SAVG  =  AbSF(PKMAG(NPK,l)  ) 

IF!  SAVc-PKMAX2)75,7(),  70 

7C  PKM AX2=SAVE 
PKTIM2=T( IMl) 

75  CONTINUE 

I F ( SENSE  LIGHT  4)77,82 
77  SPKMAX=E0(ITM1MI,1) 

PKMAX1=ABSF(SPKHAX) 

PKTIMI=T1 ITM1MI) 

NPK 1=NPK 
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78 


! 

I 


00  80  I = I TM 1 , I TM2 
I M 1  =  1- 1 


7  9 


80 


82 


PRE  V-PKC  S 

E0( I , 1 ) =E SUM < T<  I  )  ) 
PRES=EO(I*l)— EC(IM1»1) 
1F(PKEV*PRES) Tu,79,80 
NPK=NPK+1 

PKMAGCNPK,  1  )=EC(  I  M 1  ,  1  ) 
PKTIM(NPK,1)=T(  I  M 1  > 
SAVE= ABSF  ( t()  (  IM1,  1  )  ) 
IFJSAVE-PKMAX1 >80,79,79 
PKMAX 1 =S AVE 
PKf IMl=r< IM1) 

SPKMAX=E0( IH1 ,  1  ) 
CONTINUE 
I  1  =  1 TM2P1 
1 2=  I  T  1 2 
NPK2=NPK+l 
GO  TO  94 
I PERM= 1 
I TtMP=2 


RATI01=PKMAX1/PKMAX2 

AMPRNG=AMPC/PKMAX1 

IF  (RAT  10 1-10.0)  1.000*  1001,  ICO* 

1000  ARNG10=AMPRNG 
GO  TO  1002 

1001  ARNG10=10.0*AMPRNG 

1002  WRITE  OUTPUT  TAPE  6,600 

WRITE  OUTPUT  TAPE  6 , 605 , PKMAXi , PKTI Ml ,  PKMAX2 , PKT I  M2 , K AT  1 0 1 
WRITE  OUTPUT  TAPE  6 ,6 20 , ( BP (1) , I =N 1 , N2 , N1 NC  ) 

WRITE  OUTPUT  TAPE  6 , 6  1 C  ,  ( TAUPI I ) , I  =  NI»N2,N1NC) 

CALL  TICKS 
CALL  PLOT  1 

C  THE  NEXT  SECTION  PERMUTES  ALL  BIN),  ONE  AT  A  TIMt,  TESTS  FOR  ANY 
C  IMPROVEMENT,  PLOTS  IMPROVEMENT ,  PRINTS  OUT  NEW  BIN) 

84  RAT I0*UATI01 
I  CNT=  ICNT+- 1 
IFLAG-1 

DO  290  J=N1 , N2 , N I NC 
J1»J 

BTI J)-BP( JI+BINC 
IF(BINC)24C,260,250 
240  IFIBT(J) 1245,260,260 
245  BT ( J ) *0. 0 
GO  TO  260 

250  IF(BT(J)-1. 0)260, 260, 255 
255  BT ( J ) *1 . 0 

260  CALL  PKCALC(BPIJ) , rt T ( J  )  ) 

I F ( SENSE  LIGHT  1)290,265 
265  BTI J)-BPI JJ-BINC 

IFIBINC 1270,285, 290 
270  IF ( BT ( J )-l. 0)285, 285, 275 
275  8T(J)«1.0 
GO  TO  2B5 

280  IF(HTIJ) 1282,285,285 
282  BT ( J )*0.0 

285  CALL  PKCALCIBPI J) ,BT( J)  ) 

290  CONTINUE 

AMPRNG® AMPC/PKMAX 1 

IF (  RAT  101-10. 0)291, 292, 292 

291  ARNG10-AMPRNG 
GO  TU  293 

292  ARNG10=10.0*AMPKNG 

293  CALL  TICKS 
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CHNG=RAT 101-RAT  10 

1 F ( CHNG-ABSF  (  CHANGE *R A T  10  )  1300,300,295 
295  CALL  PL0T2 

WRITE  OUTPUT  TAPE  ft  ,  660  ,  PKM  AX  1 ,  PKT  l  Ml ,  PKM  AX  2 ,  PKT  1  M2 ,  k  A  I  II)  1 ,  CHNG 
_ £  WRITE  OUTPUT  TAPE  6 , 62'i ,  (  BP  II)  ,  I  =N1  ,M2  ,  N I  NC  ) 

C  THE  NEXT  SECTION  PERMUTES  ALL  TAU(N),  ONE  AT  A  TIME,  TtSTS  FOR  ANY 
C  IMPROVEMENT,  PLOTS  IMPROVEMENT ,  PRINTS  OUT  Ni-w  TAU(N) 

RAT IO  =  RA  I  lOl 
I FL AG=C 
J  1  =Nl 

T  AUT  ( J 1 ) =T AUP ( J I ) ♦TAUINC 
I F  (  PKOO  )  lift ,  1 1 0  *  1 1  3 
8ft  I  F  (  TAli  I'lC )  90  ,  llv»,  IO0 
90  IF (TAUT  I Jl)-TAuMIN)9ft,l  10,110 
95  TAUTt J1 )=TAUMIN 
GO  TO  110 

100  IFITAUT  (  J1  J-TAUMAX)  110,110,105 

105  T AUT ( J l )  =  T AUMAX 

11C  CALL  PKCALC(  TAUP( J1 ) , TAUT ( J1 I ) 

I F ( SENSE  LIGHT  3)150,115 
115  TAUT(Jl)  =TAUP  I J 1 ) -TAU INC 
IF(PR00)14ft, 145,120 
120  IF(TAUINC)125, 145,135 
125  I  F  (  TAUTT  JD-T  AUMAX)  145,145,130 
130  TAU T l J l ) =TAUMAX 
GO  TO  145 

135  I F ( TAUT ( J 1  )-T AUMIN )  140,145,145 

140  TAUT! Jl )=TAUMIN 

145  CALL  PKCALCI  T AUP( J1 ) , TAUT ( Jl ) ) 

150  00  lftO  J=N1P1,N2MI,NINC 
J1*J 

T AU f ( J ) ’TAUPI  J ) ♦ TAU I NC 
CALL  PKCALCI  T  AUPI J) , TAUT ( J ) ) 

I  F ( SENSE  LIGHT  3)160,155 
155  TAUTT J )*TAUPI  J ) -TAU INC 

CALL  PKCALCI  TAUPI  J )  ,  TAUT  ( .1 ) ) 

160  CONTINUE 
J 1  =N2 

TAUT(Jl )  =  TAUP  I Jl) +TAUINC 
IF(PROU) 190,  190,166 
165  I F ( TAUI NC ) 180 , 190 , 1  70 
170  IFITAUT (Jl)-TAUMAX) 190,190,)T5 
175  T AUT I J 1 )  =  TAUM A X 
GO  TO  190 

18U  IF ( TAUT (Jl)-T AUMIN)  185,19  0,190 

185  TAuTI Jl )=TAUMIN 

190  CALL  PKCALCI  T  AUP  I  Jl), TAU  TUI)) 

I  F  I  SENSE  LIGHT  3)230,  195 
195  TAUT(J1)*TAUPI  JD-TAUINC 
I  F ( PROO 1200,226,226 
200  IFITAUINC 1215, 225,206 
205  IF  «  TAUI  I  JD-TAUHIN)  210,226,225 
210  T AUT I J 1 ) *TAUM  1 N 
GO  TU  225 

215  IF (TAUT! Jll-TAUMAX 1225,225, 220 

220  TAUTt Jl )=TAUMAX 

225  CALL  PKCALCITAUPIJ1 ) , TAUT  I Jl ) ) 

230  AMPRNG=AMPC/PKMAX l 

IF  I  RAT  101-10.0)23 1, 232, 232 

231  ARNG10S AMPRNG 
GO  TU  233 

232  ARNG10= 10.0* AMPRNG 

233  CALL  TICKS 
CHNG=RAT I01-RAT10 

1  F ( CHNG-ABSF  I  CHANGE  *R AT 1U ) ) 300 , 300, 235 
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235  CALL  PLOT 2 

WRITE  OUTPU1  TAPE  6 , 65(1 ,  PKMAX  1 ,  PK  T I  M 1  ,  PKMAX2  ,  PK  T I  M2 ,  k  A 1  I  (H  ,  CHNG 

_ ^WKITE  OUTPUT  TAPP  6 ,  MO ,  (  T  AUP  (I)  *  I  =N1 ,  N2,  N1 NU 

l  F  <  lCNT-MAXCNDfl.4,310,310' 

300  WRITE  OUTPUT  TAPE  6 , *30 , I CNT , CHNG , I FL AG 

305  WRITE  OUTPUT  TAPC  6 , 605 , PKMAX 1 , PK T I M 1 , PKM AX2 , PK T I  M2 , RA T  1 0  1 
WRITE  OUTPUT  TAPE  6  ,62P,!  BP  II)  ,  I  =  N1 ,  N2 ,  NINC  ) 

WRITE  OUTPUT  TAPE  6 , 6 10 , < T AUP ! I) , I =N 1 , N2, N1 NC I 
CALL  PL0T1 
GO  TO  315 

310  WRITE  OUTPUI  TAPE  A , 64C ,CHNG, MAXCNT 
GO  TO  305 
315  CONTINUE 
500  FORMAT  <  7F8. 2 ) 

510  FORMAT! AI5/5E1A. A) 

520  FORMA  T 1 5E 14 . A ) 

530  FORMAT! 10F6.3) 

600  FORMAT!  A8H1INITIAL  CALCULATIONS  GIVE  THE  FOLLOWING  RESULTS/) 

605  FORMAT ! 18M0MAIN  LOBE  MAG.  =  ,  E 15 . 8 , 2X , 7HAT  T  =  , F8 . 3/23HOMAX .  SIDE 
1  LCHE  PAG.  =  ,E15.0,2X,7HAr  T  *  , F8 . 3 , 5X , 1 7HG I V 1 NG  A  RATIO  =  ,E15. 
28) 

610  F0HMAT!21H0THE  VALUES  OF  T AU ( N ) / ( 10E 1 3 . 6 ) ) 

620  FORMAT!  19H0THE  VALUES  OF  H ! N ) / ( 20F A . 3  )  ) 

A  30  FORMA  T ( 7HU I N  THE.15.30HTH  RUN,  THE  CHANGE  IN  RATIO  =  , E 1 3 . 6 , 34HW I T 
1H  THE  LAST  PER  ! URGED  SET  BEING  ,11) 

640  FORMAT! 39H0EXECUT ION  TERMINATED  WITH  A  CHANGE  OF  , E 1 3 . A ,  5 X , 6HAF T ER 

1  , I5.2X .4HRUNS ) 

650  FORMAT! IH0/18HCMAIN  LUBE  MAG.  =  , E 1 6. 8, 2X , 7HAT  T  =  ,F8.3/ 

123H  MAX.  SIDE  LUBE  MAG.  =  , E 1 5 . 8 , 2X , 7H A T  T  =  .F8.3/13H  NEW  RATIO  = 

2  , l l 3 . 6 , 5 X , 1 2h A  CHANGE  UF  ,EI3.6) 

CALL  EX  I  I 

ENO!  1,0, u, 0,0,0, l ,0,0,0,0,0,0,0,01 


4 .  Subprogram  ESUMI 


The  following  function  is  written  as  a  subprogram  to  allow  changes 
in  the  scan  rate  and  the  initial  frequency  to  be  made  without  the  need 
to  recompile  the  main  program: 

m 

ESUMI  =  Y  A(T-T  )  B  cos  <  2n  2 (T-T  )  -  0.0l(T-T  )2 
Z-  n  n  \  L  n  n 

n=l 


where  T  is  the  current  value  of  t  in  the  main  program. 


t 

4 

f 

B- 


u 


ESUM  SUBPROGRAM  TO  CALCULATE  SUM  OF  COSINES 
FUNCTION  ESUM! TIME) 

DIMENSION  T( 2501  ),TAUP( 101 J.TAUT!  101 ) ,BPI 101 ) ,BT 1 101 ) ,E0!2501 , 2 ) « 
1PKMAGI2000, 2 ) , PKTIM! 2000 ,2  I ,PKRAT(  2000 ) , AT 1CK ( 20) , GT ICK (20 ) 

COMMON  TMAXI , T  INC ,  TR, TLNGTH.TRNGE , TLAG, LASTT, IT  12, I TM1 , 1TM2  , 
1ITM1M1, ITM2PI,NI,N2,NINC,NPK,NPK1 ,NPK2 , PKMAXl , PKMAX2.PKTIM1 , 
2PKTIM2.IPERM, I  TEMP , J 1 , 1 FL AG, AMPRNG, A 1 , A2, ARNG10.RAT 1 01 , SPKMAX, 

3T  ZERO, T  ,TAUP,TAUT,HP,6T,E0*PKMAG»PKT1M»PKRAT,ATICK»BTICK»TL»TL 1 
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ESUM-0.0 

00  10  J»N1,N2,MNC 
THTAU*T  I  HE-TAUPI J I 

ESUM»ESUM*BP«  J)*AC  TMAU)»C0SF  <6. 2831853*  12. 0*TMTAU-0.01«TMTAU* 
lTMTAUn* 

10  CONTINUE 
RETURN 

END! 1,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0) 


5 .  Subprogram  ESUMII 

ESUMII  differs  from  ESUMI  by  the  addition  of  a  phase  term  to  the 
cosine  argument.  The  function  calculated  is 


m 

ESUMII  =  ^  A(t-T^)  Bn  cos  |2n|^2(T-Tn)  -  0.0l(T-Tn)2  +  PHII 
n  ^1 


where  T  is  the  current  value  of  t  in  the  main  program  and  PHII  is 
a  phase  term  obtained  from  the  main  program.  ESUMII  can  only  be  used 
with  Program  I . 


ESUM  SUBROUTINE  TO  CALCULATE  INITIAL  FUNCTION 
FUNCTION  ESUMI TIME! 

DIMENSION  T  (  2*501  >,  TAUP  (  101),  TAUT!  101),  BPI  10 1 ) ,  BT 1 10 1 ) ,  E0 1  250 1 , 2  )  , 
IPKMAGt 2000, 2) ,PKTIM( 2000,2 ),PKRATI 2000) .AT1CKI 20), BTICKI20), 
2PHII21) 

COMMON  TMAXI , TINC.TR, TLNGTH, T RNGE , TL AO , LAST T , I T 1 2 , ITM1, ITM2, 
1ITM1M1,  I  TM2P1 , Nl , N2 ,N INC « NPK, NPK1  ,NPK2 , PKMAXl , PKMAX2 , PKT IM1 , 

2PK1 IM2, IPERM, ITEMP, Jl,  l  FL AG, AMPRNG,  A 1 , A 2, ARNG10 ,R A TIOl.SPKMAX, 

3  T  ZERO, T  ,  T  AUP , T AUT , BP , 81 , EO , PKMAG , PK 1 1  M , PKR A T , AT  I CK , B T I CK , TL , TL 1 , 
4PHI.PHII 
ESUK=0.0 

00  10  J=N 1  *  N2 » N INC 
TMTAU=I 1MC-TAUPI J) 

01  *BPI J 1 • A  C TMTAU) 

B2=2. 0«TMT  All-0. 0 1  »TM1  AU»TMTAU 
ESUM*ESUM«-rtl*CUSFC6.2B3lB,5  3»IB2  +  PHI  I  ) ) 

10  CONTINUE 
RETURN 

END! 1, 0,0, 0,0,0, 1 ,0,0, 0,0, 0,0, 0,0) 


69 


SEL-65-043 


6 .  Subprogram  EREV 

After  a  perturbation  of  B  or  T  has  been  made,  EREV  calculates 
the  new  value  of  r(t)  for  all  values  of  t  in  the  main  program.  The 
new  value  is 


.  /  1\  1 

L 

„  /  i\ 

/  !\2 

1 

+  A (t-Tj  cos  j 

[2lT 

-  O.Ol(t-T^) 

1 

-  A(t-1\)  B^  cos  | 

[2” 

2(t-T.) 

-  0.0l(t-Ti)2 

!} 

where  r  (t)  is  the  value  of 
o 


S  rn(t)  =  r(t) 


i=l 


before  the  perturbation.  The  change  in  r(t)  due  to  the  change  in  the 
th 

i  tap  is  calculated  by  subtracting  the  initial  contributions  from  the 

i  tap  from  rQ(t)  and  adding  the  new  contributions.  This  must  be  done 

for  each  value  of  t . 


C  EREV  SUBROUTINE  TO  REVISE  FUNCTION 
SUBROUTINE  EREV 

DIMENSION  T(2SU1)*TAUP(101)#TAUT(101)*BP(101)*BT(101)*EO(2501*2)* 
lPKMAG(2000»2),PKTIM<2000*2)»P<RAT(2000>»ATICK(20)*eTICK(2o> 

COMMON  TMAXI*TINC*TR»TINGTH*TRNGE»TLA(5,1*STT*ITI2#ITMI,ITM2, 
IITMIM1# ITM2P1,N1#N?#N1NC*NPK*NR«!,NPK2,PKM*X1»PKM»X2#PKTIM1# 

2PKTIM2#  I  PERM,  I  TEMP#  J1*IFLAG*  ANPRNG, Al, A2»ARNGi0*RATI0l»SPKMAX, 
3TfER0»T#TAUP* TAUT* 8P*BT»E0»PKMAG»P«TIM»PKRAT»ATIC<»BTICK»TL» Til 
Bl-BPCJl) 

TAUI-TAUPU1) 

I F  C I  FLAG  )5*  5* 1 0 
5  B2>B1 

TAU2»  T  A  U  T ( J 1 ) 

GU  TO  IS 
10  82»BT ( J 1 ) 

TAU2-TAU1 
15  00  20  II«1*ITI2 
TMTAU1»T(II)-TAUI 
TMTAU2»TCII)-TAU2 

E0(II#ITEMP)«E0(II*IPERM)-BI*A(TMTAUI)*C0SF{6.2B31653*(2.0*TMTAU1- 

l0.0100*TMTAUi*TMTAUl))t82*A(TMTAU2)*COSF(6,283l853*(2.0*TMTAU2- 

20.0100*TMTAU2*TMTAU2)) 

20  CONTINUE 
RETURN 
END 
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7 .  Subprogram  TAUN 


If  ITAU  =  0,  then  this  subprogram  is  called.  It  calculates  the 
values  of  T  which  will  cause  the  contributions  from  all  taps  to  add 
in  phase  at  t  =  50  sec  as  described  in  Eq.  (5).  The  function  calcu 
lated  is 

TAU  =  100. 0\/ 1.01  -  (O.Ol)FLN  -  50.0 

where  FLN  is  a  floating  point  N. 

TAUN  SUBPROGRAM  TO  CALCULATE  TAU  AS  A  FUNCTION  OF  N 
FUNCTION  TAU I FLN I 

TAU*100.0*SQRTF.(  1.0l-0.01»FLN)-50.0 
RETURN 

ENDI 1 ,0, 0, 0. 0. 0, 1. 0,0,0, 0, 0,0, 0.0) 


8.  Subprogram  AN 

This  subprogram  calculates  the  values  of  A(t-Tn). 


t-T  -25\ 

— 55 — )  -io  *  (t-Tn)  *  60 


0 


-10  >  (t-T  )  >  60 
n 


AN  SUBPROGRAM  TO  CALCULATE  A  AS  A  FUNCTION  OF  T-TAUIN) 

FUNCTION  A(TMTAU) 
lFITMTAU+10.0) 15,10,5 
5  IFITMTAU-60.01 10,10, 15 
10  A«1.0-.81632653E-3*ITMTAU-25.0)*ITMTAU-25.0) 

RETURN 
15  A«0.0 
RETURN 

EN01 1,0, 0,0, 0,0, 1,0,0, 0,0, 0,0, 0,0) 
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9 .  Subprogram  PL0T1 


PLOT1  is  a  subprogram  which  plots  r(t)  using  Subprograms  XYPLOT, 
,<YLABL,  and  XYFRC.  The  largest  peak  of  r(t)  is  scaled  to  4.5  in. 
which  places  r(t)  =0  along  the  center  of  the  plot.  PLOT1  calls 
PL0T2  before  returning  to  the  main  program. 


PL0T1  SUBPROGRAM  TO  PLOT  INITIAL  ANO  FINAL  DATA 
SUBROUTINE  PL0T1 

DIMENSION  T  (2  501 ), TAUP! 101 J,  TAUT!  1 01 ) ,BP! 101 ) , BT ! 101 )  , 60 ! 2 501 , 2 ) , 
1PKMAG<2000,2) , PKT l M! 2000,2 ),PKR AT! 2000) , AT ICK!  20 ) , BT ICK I  20 ) 

COMMON  TMAXI , T I NC , TR, TLNGTH, TRNGE .TLAG.LASTT, IT  12, ITM1 , 1 TM2 , 
1ITM1M1, I TM2P1 ,N1,N2,N1NC,NPK,NPK1 ,NPK2, PKMAX1 , PKMAX2 , PKT I M 1 , 
2PKTIM2, IPERM.ITEMP.Jl , I  FLAG, AMPRNG , AI,A2,ARNG10, RAT  101, SPKMAX, 
3TZER0.T.TAUP, TAUT, BP.BT ,E0 .PKMAG , PKT l M, PKRAT , AT  I CK ,BTI  CK, TL.TLl 
CALL  AXP LOT (AMPRNG ) 

CALL  XYPLOT IT ( 1 » ,E0< 1, I  PERM)  ,  TRNGE , AMPRNG, TZE RO, 0 . 5 , TLNGTH , 3) 

DO  20  1*2, I TI 2 

CALL  XYPLOT! T< I  1 ,E0! I , I  PERM )  ,  TRNGE , AMPRNG, TZE RO, 0. 5 , TLNGTH, 2 ) 

20  CONTINUE 
140  CALL  XYFRC 

CALL  XYPLOT! 1 .0,0.0, 1.0,1.0,0.0,0.0,1.0,31 
CALL  XYFRC 
CALL  PLOT2 
RETURN 

END!  1,0, 0,0, 0,0, 1,0, 0,0, 0,0,0, 0,0) 


10.  Subprogram  PL0T2 

PLOT2  is  the  same  as  PLOT1  except  that  the  largest  peak  is  scaled  to 
45  in.  instead  of  4.5  in.,  and  the  region  from  1MAIN1  to  TMAIN2  is  not 
plotted . 


PL0T2  SUBPROGRAM  TO  PLOT  INTERMEDIATE  DATA 
SUBROUTINE  PL0T2 

DIMENSION  T ( 2501 ) , TAUP (101)»TAUT!101)»HP! 101 ) ,bT 1 1011 ,  EO!  2501 , 21 , 
1 PKMAG ! 2 000,2) , PKT IM 12000,2 ) , PKRAT ! 2  000  I , AT ICK ( 20 1 , BT ICK ! 20 ) 

COMMON  TMAXI ,TINC, TR, TLNGTH, TRNGE, TLAG, LASTT, I TI2.ITM1, I TM2, 
1ITM1M1, ITM2P1,NI,N2,NINC,NPK,NPK1 , NPK2 , PKMAX 1 , PKMAX2, PK T I M 1 , 
2PKTIM2, 1  PERM, I  TEMP, Jl, I  FLAG, AMPRNG, Ai,A2,ARNGlO, RAT  101,  SPKMAX, 
3TZER0, T, TAUP, TAUT, BP, BT,E0,PK MAG, PK TIM, PKRAT , AT  I CK, BT ICK , TL , TL 1 
CALL  AXPLOT ( ARNG10 ) 

CALL  XYPLOT! T! 1 ) , EO! 1 , 1  PERM) ,  TRNGE , ARNG 10, TZERO, 0. 5, TLNGTH, 3 ) 

DO  10  L* 1 , NPK 1 

CALL  XYPLOT  (PKTIM!L, I  PERM) , PKMAG < L , I  PERM) , TRNGE , ARNG 10,  TZERO, 0 .5, 
1TLNGTH , 2  ) 

10  CONTINUE 

ET1*E0< ITM1MI, IPERM) 

ET2=E0I ITM2P1, IPFRM' 

E  T  3*ET l * ARNG1 0  ♦ 

ET4=ET2*ARNG10«-0.5 

WRITE  OUTPUT  TAPE  6 , 2000, PKTI  M ( NPK l , I Pt RM) , PKMAG! NPKl , IPERM ) , 
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IT  (  ITM1M1 )  ,£Tl,ET3,  TUTM2P1)  ,E  T2  ,E  T4  ,PKUM(  NPK2, 1  PERM) , 

2PKMAGI  NPK2 »  I  PERM) 

2000  FORM*  T ( 10E12.4) 

IF(ET3-1. 0)12. 12.14 

12  CALL  XYPLOTm  I TM1M1 > , ET l , TRNGE , ARNG10, TZE R0.0.5, TLNGTH, 2 ) 

14  IF(ET4-1.0)16, 16, 18  „ 

16  CALL  XYPLOT(T< 1 TM2P1 ), ET2 , TRNGE .ARNGIO, TZERO, 0.5 , TLNGTH,  3 ) 

GO  TO  19 

18  CALL  XYPL0TJPKTIM1NPK2.IPERM) , PKMAG ( NPK2, 1  PERM ) , TRNGE,  ARNG10, 
1TZER0.0.5, TLNGTH, 3 ) 

19  DO  20  L«NPK2,NPK 

CALL  XYPL0T1PKTIMIL.IPERM) , PKMAG ( L. I  PERM) , TRNGE, ARNGIO,  TZER0.0.5, 
1TLNGTH.2 ) 

20  CONTINUE 

CALL,  XYPLOTm  1  T 1  2  ) ,  EOI ITI 2 , 1  PERM) ,  TRNGE  ,  ARNG  10,  TZERO.O.  5, 
1TLNGTH.2) 

CALL  XYFRC 

CALL  XYPLOT<1.0,0.0,1.0,1.0,0.0,0.0, 1.0,3) 

CALL  XYFRC 
RETURN 

END( 1,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0) 


11 .  Subprogram  PKCALC 


This  subprogram  evaluates  the  new  set  of  maxima  and  minima  after  a 
perturbation  has  been  made.  If  an  improvement  in  the  peak-to-sidelobe 
ratio  has  been  made,  PKCALC  returns  to  the  main  program  and  leaves  the 
variable  changed.  If  no  improvement  is  found,  the  sign  of  the  perturba 
tion  is  changed  and  the  peak-to-sidelobe  ratio  is  again  calculated.  As 
before,  an  improvement  will  return  control  to  the  main  program;  if  no 
improvement  is  obtained,  the  variable  is  returned  to  its  original  value 
before  control  is  returned  to  the  main  program. 


PKCALC  SUBPROGRAM  TO  CALCULATE  NEW  SET  OF  PEAKS  AND  EVALUATE  THEM 


SUBROUTINE  PKCALC 1 AP, AT ) 

DIMENSION  T( 250 1 ) , TAUPC 101 ) ,TAUT( 101) ,BP( 101 ) ,BT ( 101 ) ,E0(2601 ,2 ) » 
1  PKMAG  1 2000, 2)  *  PKT I Ml  2000,2 )  ,PKRAT(2000 )  ,AT  ICK 1 20) ,  DTICM20) 

COMMON  TMAXI.TINC, TR, TLNGTH, TRNGE , TLAG, LAST! ,  ITI2, ITM1, ITM2, 
1ITM1M1,  ITM2P1 ,N1,N2,NINC,NPK,NPK1,NPK2, PKMAXi  , PKM AX2 , PKT IM1 , 
2PKTIM2, I  PERM, I  TEMP , J1 .IFLAG.AMPRNG, Al, A2, ARNG10.RAT 101, SPKKAX, 
3TZER0, T, TAUP, TAUT, bP, BT,E0, PKMAG, PKT 1M.PKRAT,  ATI CK.BT ICK, TL.TLl 
I F ( SENSE  LIGHT  3)5,5 

5  NOPKS-O 
CALL  EREV 

PRES«E0(2, ITEMP)-E0( 1,1  TEMP) 

1F(A8SF(EO(2,IT£MP))-ABSF(EO(1,I TEMP) ) ) 7,6,6 

6  TPMAX2*ABSF (EO 1 2, 1  TEMP ) ) 

TPTIM2*TI2) 

GO  TO  8 

7  TPMAX2* ABSF (E0{ 1,1 TEKP) ) 

TPTIM2*T (  1 ) 
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6  SENSE  LIGHT  4 

no 

I  2= I TM1M 1 
9  00  29  1*11,12 
I  Ml* I- 1 
PREV*PRES 

PRES-EO!  I,lTEMPt-EO(IMl,  I  TEMP) 
IF(PREV»PRES> 10,10,29 
10  NOPKS-NOPKS+1 

PKMAG (NOPKS, I  TEMPI »E0<  IMl.lTEMP) 
PK TIM! NOPKS  » ITEMPI*T ( IM1 ) 

SAVE* ABSF ( PKMAG INOPKS,  l TEMP  H 
I F (SAVE- TPM AX 2 129,25,25 
25  TPMAX2-SAVE 
TPTIM2*T ( IM1) 

29  CONTINUE 

IF (SENSE  LIGHT  4130, 34 

30  STPMAX*EO( I TM1M1, I  TEMP  I 
TPMAX1*ABSF ( ST PM AX) 

TPTIM1-T! ITM1M1I 
NPKT1«N0PKS 

00  33  1*ITM1,  ITM2 

IM1-I-1 

PREV*PRES 

PRES*EO ( I , I  TEMP l-EOI  INI ,  I  TEMP  1 
IF(PREV«PRES)31,31,33 

31  NOPKS-NOPKS+1 

PKMAG (NOPKS, I  TEMP l*EO(  I  Ml, I  TEMP  I 
PKT IM( NOPKS, ITEMP|*T(IM1I 
S  AVE*AbSF ( PKMAG! NOPKS, 1  TEMP  1 ) 

I F ( SAVE -TPM AX  1 133,32,32 

32  TPMAX1*SAVE 

ST PM AX* PKMAG! NCPKS.ITEMP) 
TPTIM1*T( IM1I 

33  CONTINUE 
I1*ITM2P1 
12*  IT  12 
NPKT2-N0PKS+1 
GO  TO  9 

34  TRAT10*TPMAX1/TPMAX2 

IF (TRA1I0-RATI01 150, 50,35 

35  PKMAX1*TPMAX1 
PKTIM1*TPTIM1 
PKMAX2*TPMAX2 
PKT IM2-TPT IM2 
RATI01«TRATI0 
SPKMAX*STPMAX 
NPK-NOPKS 
NPK1 *NPKT1 
NPK2*NPKT2 
AP«AT 

SENSE  LIGHT  3 
I FIIPERM- 1140,40,45 
40  IPERM-2 
ITEMP-l 
GO  TO  50 
45  IPERM-1 
I  TEMP-2 
50  CONTINUE 
RETURN 

END! 1,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0 


SEL-65-043 


74 


12.  Subprogram  AXPLOT 


This  subprogram  draws,  a  box  around  the  plot.  AXPLOT  is  called  by 
PL0T1  or  PL0T2 ,  and  AXPLOT  calls  TICKS. 


AXPLOT  SUBPROGRAM  TO  PLOT  AXES  FOR  GRAPHS 
SUBROUTINE  AXPLOT ( ARNG ) 

DIMENSION  T( 2501 ) ,TAUP 1101 > ,TAUT( 101), HP! 101) ,BT( 101 ) , EOI 2501 • 2 ) , 
1 PKMAGI 2000 ,2 ) .PKTIM1 2000,2 ) ,PKR AT  12000)  .ATICKI 20 ) , BT I CK (20 ) 

COMMON  TMAX1 ,T INC, TK,TLNGTH,TRNGE ,TLAG, LASTT, IT  12, I TM1 , ITM2, 
1ITM1M1,ITM2P1.N1,N2,NINC,NPK,NP>U,NPK2,PKMAXI,P<MAX2,PKTIM1, 
2PKTIM2,IPERM,ITEMP,J1, IFLAG, AMPRNG, A1 , A2, ARNG10,  RAT IQl.SPKMAXf 
3TZER0,T,TAUP,TAUT,BP,BT,E0,PKMAG,PKTIM,PKRAT,ATICK,BTICK,TL,TL1 
CALL  XYPLOTI 1.0, 0.0, 1.0, 1.0, 0.0, 0.0, 1.0, 3) 

CALL  XYFRC 

CALL  XYPLOT(0.0,0.0,TRNGE,0.1 ,0.0,0.5, TLNGTH.3) 

Tl-TLAG-10.0 
00  10  L-l, LASTT 
T1-T1M0.0 

CALL  XYPLOT ( T 1 , 0. 0 , TRNGE ,0. 1,0.0 , 0. 5, TLNGTH, 2 ) 

CALL  XYPLOT (T1 , 0.05, TRNGE ,0. 1 ,0.0, 0.5, TLNGTH, 2 ) 

CALL  XYPLOT <T1 ,-0.05, TRNGE ,0. 1 ,0.0,0. 5, TLNGTH,  2) 

CALL  XYPLOT 1 T 1,0.0, TRNGE, 0. 1,0. 0,0. 5, TLNGTH, 2) 

10  CONTINUE 

CALL  XYPLOT ITR* P.0, TRNGE, 0.1, 0.0, 0.5, TLNGTH, 2) 

CALL  XYPLOTIO.O.Al.l. 0,0.1 ,0.0, 0.0, TLNGTH, 3) 

00  15  L-1,10 

CALL  XYPLOTI 0.0, ATICKI L) , TL , ARNG.O. 0, 0 . 5, TLNGTH, 2 ) 

CALL  XYPLOTI-0.05.ATICKIL) ,TL, ARNG, 0.0,0. 5, TLNGTH, 2 ) 

CALL  XYPLOTI 0.0, ATICKI L ) ,TL. ARNG, 0.0,0. 5, TLNGTH, 2) 

15  CONTINUE 

CALL  XYPLOT I 0.0.A2, 1.0,0. 1 ,0.0, 0.0, TLNGTH, 2) 

CALL  XYPLOTI 1.0, A2. 1.0, 0.1 ,0.0, 0.0, TLNGTH, 2) 

00  20  L-1,10 

CALL  XYPLOTI TLNGTH, BTICKI L ) ,TL, ARNG, 0.0, 0.5, TLNGTH, 2 ) 

CALL  XYPLOTI TL1, BTICK IL ) ,TL, ARNG, 0.0, 0.5, TLNGTH, 2) 

CALL  XYPLOTI  TLNGTH.BTICML  )  ,TL ,  ARNG,  0.0 ,0. 5,  TLNGTH,  2 ) 

20  CONTINUE 

CALL  XYPLOT 1 1.0, A  1,1. 0,0. 1,0. 0,0.0, TLNGTH, 2) 

CALL  XYPLOT I  0.0, A  1,1. 0,0.1 ,0.0, 0.0, TLNGTH, 2) 

RETURN 

ENOI 1,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0) 
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13.  Subprogram  TICKS 

TICKS  draws  tick  marks  on  the  side  of  the  box  drawn  by  AXPLOT  to 
mark  0.1,  0.05,  0.025,  0.0125,  and  0.01  times  the  maximum  value  of 
r(t). 

TICKS  SUBROUTINE  TO  CALCULATE  AMPLITUDE  TICK  MARKS 
SUBROUTINE  TICKS 

DIMENSION  T(2501),TAUP(101 I, TAUT! 1011, BP! 1011 ,BT ( 101 ) , E0( 2501 , 2 ) , 

1 PKMAGI2000. 2) ,  PKTIMt  2000*2  t.PKRATI  20001  •ATICK1201  , BTICKI20) 

COMMON  TMAXI , T INC , TR, TLNGTH.TRNGE , TLAG, LASTT , ITI2, ITM1.ITM2, 
lITMlMl(ITM2Pl*NltN2*NINC(NPK*NPKl«  NPK2 • PKMAX 1  * PKMAX2*  PKT IM1 * 

2 PKT I M2* I  PERM* I  TEMP  * J1 * 1FLAG* AMPRNG* A1 * A2*  ARNG10* RATIOltSPKMAX* 
3TZER0* T, TAUP. TAUT, BP tBT.EO.PKMAG, PKT IM.PKRAT, ATICK, STICK, TL.TL1 


ATICKI 10 

) *0. 1 »PKMAX1 

ATICKI9) 

■0 . 05  »PKMAX 1 

ATICK(B) 

«0.025»PKMAX1 

ATICKI7) 

■0.0125»PKMAX1 

ATICKI6) 

■0. 01 *PKHAX1 

ATICKI5) 

— ATICKI6) 

ATICKI4) 

— ATICKI7) 

ATICKI3) 

— ATICKI8) 

ATICKI2) 

— ATICKI9I 

ATICK(l) 

■-ATICKI 10) 

BTICK(l) 

■ATICKI 101 

BT ICK ( 2 ) 

■ATICKC9) 

BTICKI3) 

■ATICKI 8) 

BTICKI4) 

■ATICKI7) 

BT ICK ( 51 

■ATICKI6) 

BT ICK ( 6 ) 

■ATICKI5) 

BT ICK ( 7 1 

■ATICKI 4) 

BT ICK (81 

■ATICKI3) 

BTICKI9) 

-ATICKI2) 

8TICKI10I-ATICKII) 

RETURN 

END! 1*0* 

0*0*0* 0*1 *0,0 *0*0 *0*0*0,01 
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